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DESIGN AND PERFORMANCE AT A LOCAL MACH NUMBER OF 6 
OF AN INLET FOR AN INTEGRATED SCRAMJET CONCEPT 

Carl A. Trexler and Sue W. Souders 
Lai^ley Research Center 

SUMMARY 

A research program on hypersonic propulsion at the NASA Langley Research Center 
is focused on the development of a concept for a modular supersonic combustion ramjet 
(Langley Seramjet Module). The modular engine concept is designed to integrate with the 
airframe; precompression of the ermine airflow will be produced by the vehicle bow shock 
and additional expansion of the nozzle exhaust will be produced by the vehicle afterbody. 

As part of this research program, component Investigations are in progress on the base- 
line inlet configuration and the present paper reports the design philosophy and results of 
experiments at Mach 6 to evaluate the performance of the inlet. 

With the integration advantages, the inlet was designed with modest contraction 
ratios and fixed geometry. Three fuel injection struts contribute to the inlet flow com- 
pression and provide a short combustor design that results in low internal cooling require- 
ments. The baseline inlet configuration is rectangular in cross-sectional shape, has 
sweptback sidewall planar compression surfaces, has an opening upstream of the cowl 
leading edge through which spillage occurs for startii^ and normal operation, and has the 
external cowl surface alined with the local flow to provide a minimal external drag. The 
inlet model had a projected geometric capture area measuring 19.05 cm high by 15.24 cm 
wide. The sidewalls and struts had 48° swept leadii^ edges and nominal compression sur- 
face angles of 6°. 

Tests were conducted in the Langley 20-inch Mach 6 tunnel and the inlet model was 
instrumented to obtain both wall and survey pressure measurements which were used in 
computing performance and capture flow. The data-reduction system provided integrated 
performance data as well as contour maps of parameters such as total-pressure recovery 
and Mach number in the inlet throat. The difficulty and importance of properly position- 
ing the shock waves in the throats of hypersonic inlets were demonstrated, but no adverse 
effects were noted as a result of the inlet ingesting a boundary layer on the top surface 
which simulated the boundary layer that would be ii^ested from the vehicle forebody. 

The average throat Mach number was 3.1 compared with the predicted value of 3.4. 
The kinetic energy efficiency was 97,7 percent (0.59 recovery) compared with prediction 



of 98.3 percent (0.67 recovery), which did not account for all sources of total-pressure 
loss. The average inlet aerodynamic contraction ratio was 7.0, which does not include 
the compression expected from the vehicle bow shock. The inlet captured 94 percent of 
the flow at its face; and, overall, the inlet performance was well within the acceptable 
range for high engine performance. 


INTRODUCTION 

The attractive potential of hypersonic flight with air-breathing propulsion has been 
recognized for the past 15 years; however, major advances in technology are required. 
Exploratory research on concepts for hypersonic air-breathing engines has been pursued 
in substantial research and development programs and a broad technology base has been 
established. See, for example, references 1 to 4. The investigation of several small- 
scale, hydrogen fueled, supersonic combustion ramjet (scramjet) ermine designs has 
shown that the scramjet is a feasible engine concept and practical levels of thrust have 
been demonstrated (ref. 5). Hypersonic research and technology programs have been 
conducted (refs. 6 and 7) which illustrate the next logical step in scramjet evolution, 
which is the development of engine concepts which will integrate with the airframe. Inte- 
gration includes the use of the vehicle forebody to precompress the engine airflow before 
it enters the inlet and the use of the vehicle afterbody for additional expansion and thrust 
vectoring of the nozzle exhaust gas. Other principal design criteria for hypersonic sys- 
tems are minimum engine cooling requirements to make part of the heat sink of the hydro- 
gen fuel available for active cooling of the airframe, fixed geometry to reduce weight and 
system complexity, and minimum external drag. 

Detailed analytical and experimental studies at the Langley Research Center have 
resulted in the definition of the Langley Scramjet Module, with design features in both the 
inlet and combustor which will satisfy the engine design criteria, when the benefits of 
vehicle and propulsion system integration are included. This report deals primarily with 
the design and performance evaluation for the hypersonic inlet concept for the Langley 
Scramjet Module. The design criteria were met with the use of swept compression sur- 
faces, which produced oblique shock waves, and a matching swept throat and combustor. 
Because of the complexity of the flow, it was necessary to optimize the selected inlet con- 
f^uration from the experimental results of several earlier configuration studies. Once 
the concept was derived, computer programs aided in locating shock waves within the 
inlet, determined boundary -layer eorreetions to the interior walls, provided theoretical 
performance results, and reduced experimental test data. 

A model of the inlet portion of the Langley Scramjet Module, 19.05 cm high by 
15.24 cm wide, was tested in the Langley 20-inch Mach 6 wind tunnel. This test 
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condition represents local inlet face conditions for a flight Mach number of approxi- 
mately 7.6 after compression from the vehicle forebody. The tunnel free -stream 
total temperature and pressure were 467 K and 11.9 atm (1 atm « 101.3 KN/m^), 
respectively, and provided a Reynolds number per meter of approximately 9.8 x 10®. 
One run was made at a reduced Reynolds number per meter of 3.3 x 10®. 

SYMBOLS 

When two symbols are given for the same concept, the second one is that used for 
the computer data. 

cross-sectional areapf a stream tube 
cross-sectional area of a stream tube with sonic velocity 

Al Pf A* 

aerodynamic contraction ratio, ~ 

Al pt^i A 

distance from cowl tip (fig. 13(d)) 
distance from cowl leading edge (fig. 13(d)) 
inlet height, 19.05 cm (7.5 in.) 

Mach number 
static pressure 

static pressure in front of inlet 
total pressure 
PpitotjPITOT pitot pressure 
R Reynolds number 

s distance from foreplate leading edge (fig. 13(a)) 


A 

A* 

A/Ai 

c,C 

c’,C’ 

H 

M 

P,P 

PpPl 

Pt 


3 



s’ distance from sidewall leading edge (fig. 13(c)) 

s" distance from strut leading edge (f^s. 13(e) to 13(g)) 

T temperature 

Taw adiabatic wall temperature 

Tt total temperature 

V velocity vector 

V velocity 

W throat gap or width of capture measurement station, cm 

(figs. 11 and 25) 

X distance downstream of intersection of sidewall with foreplate 

(fig. 13(b)) 

X axis parallel to free-stream flow (fig. 57) 

V distance from foreplate (fig. 13(c)) 

y axis perpendicular to free-stream flow (fig. 57) 

Z distance away from model center line (figs. 13(a) and 13(d)) 

Z' distance across throat or across duct (fig. 11) 

z axis perpendicular to free-stream flow and y-axis 

(figs. 4 and 57) 

6 boundary-layer thickness 

6n flow turning angle normal to leading edge of swept wedge (fig. 57) 

6jjy cross -flow angle (fig. 57) 

6xz wedge angle or flow turning angle in xz- plane (fig. 57) 
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6 jj 2 shock -wave angle in xz-plane (fig. 57) 

kinetic energy efficiency 

A sweep angle measured in xy-plane, deg 

X distance around capture measurement station (fig. 62) 

p density 

<p dihedral ar^le measured perpendicular to ridge line, deg 

yp dihedral parameter 

Subscripts: 

1 conditions at inlet face or ahead of a shock wave 

2,3,4 conditions behind 1st, 2d, and 3d shock waves, respectively 

n normal to leading edge (fig. 57) 

t tangential to leading e<%e (fig. 57) 

°o free stream (in front of vehicle) 

INLET DESIGN CRITERIA AND CONCEPT 
Airframe -Engine Integration 

The advantages of integrating the airframe and engine of a hypersonic vehicle are 
well known. The parametric analysis of reference 8 indicates that the contributions of 
the vehicle forebody and afterbody are responsible for up to 70 percent of the thrust. At 
hypersonic speeds very large eigine airflows are required for adequate thrust in spite of 
the high-energy potential of the hydrogen fuel. These engine airflow requirements are 
best met by utilizing the precompression obtained from the vehicle forebody and locatiig 
the engine on the underside of the vehicle toward the aft end. (See f%. 1.) A method of 
designiig the forebody is discussed in reference 6. The propulsion system inlet area is 
therefore restricted to the space between the vehicle undersurface and the bow shock and 
is several times wider than it is high, as shown in the cross section of figure 1. This 
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geometry concept suggests the arrai^ement of a number of adjacent rectangular engine 
modules, and the use of modules also permits development in ground test faeilities of 
reasonable size. The relatively thick turbulent boundary layer generated on the vehicle 
forebody is an unfavorable characteristic which must be considered in the inlet design. 

The design of the vehicle afterbody, discussed in reference 7, is also important because 
of the large gross -thrust and moment forces involved, which can generate large trim drag 
penalties if not correctly considered. 

The Scramjet Ermine 

The primary objectives of the scramjet ermine des^n are: to provide a h^h level 
of thrust and specific impulse with efficient capture over the fl^ht Mach number range 
from 3.5 to 10, to Imve low cooling requirements in order to make a portion of the fuel 
heat sink available for active cooling of the vehicle structure at high Mach numbers, to 
have satisfactory operatiig characteristics over the Mach number range including the 
establishment of supersonic flow (starting) within the inlet at the low end of the Mach num- 
ber raige, to have fixed geometry in order to reduce system complexity as well as joint 
and seal problems, to iigest successfully the vehicle -forebody boundary layer, and to pro- 
duce low external drag. Many of these objectives are interrelated and trade-offs dis- 
cussed in reference 7 indicate that a fixed geometry inlet with moderate contraction for 
starting at a low Mach number is desirable. A moderate contraction and low internal 
pressure will not only mean a reduction in engine weght and cooling requirements but 
also an increased ability for the engine to ingest the vehicle -forebody boundary layer. 

Because the use of only fuel injection from the sidewalls would produce very long 
mixing lengths for this type of modular design, the scramjet engine concept of fgure 2 
has three struts to provide six planes of instream fuel injection. This feature not only 
reduces cooliig and shortens the combustor but also the inlet, since the struts provide a 
significant part of the inlet flow compression. The sidewalls are the main inlet com- 
pression surfaces, whereas the top surface partially eliminates expansions produced by 
the downflow created by the swept shock-wave system, which is a unique characteristic 
of the swept inlet design. The cowl is kept nearly parallel to the vehicle underbody to 
minimize external drag. 


The Inlet Concept 

The inlet must efficiently compress the airflow captured for the combustor. Refer- 
ence 7 indicates that a contraction ratio between 6 and 10 would be satisfactory for a fixed 
geometry inlet at a flight Mach number of about 1 ^, and several inlet conf^rations were 
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investigated which would fit into the area provided beneath the vehicle. Inlet interior 
walls consisting of swept planar surfaces were assumed; this assumption simplifies the 
analysis and avoids the need for three-dimensional characteristic computer programs. 
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Compression angles between 6° and 8° were employed; these angles are a compromise 
between high angles which contribute to high total-pressure losses from shock waves and 
increase the possibility of shock-induced boundary -layer separation and low compression 
angles which make the inlet long. 

Prior to the development of the inlet shown in figure 2, a design utilizirg swept com- 
pression surfaces in which the top surface was the primary compression surface was con- 
sidered and is discussed in reference 9. Disadvantages to this design were: a corner 
flow problem originating at the top surface and covering much of the inlet throat; diffi- 
culty in obtaining good capture eharacteristies over the Mach number range; and no effec- 
tive way of dealing with the vehicle -forebody boundary layer. Therefore sidewalls with 
swept leadii^ edges were made the primary compression surface for the inlet in figure 2. 
Because planes of constant flow properties tend to be parallel to the sweep lines, the fuel 
injection struts and all downstream stations are also swept at the same angle. This design 
generates a system of swept shock waves which turns the flow away from the top surface 
and thus reduces the corner flow and boundary -layer problems on that surface. A cowl 
design which would provide 100-percent captvire at high Mach numbers was tested on a 
preliminary design, and those results indicated that better startirg and operating perfor- 
mance could be obtained with the pointed cowl leading edge located near the struts as 
shown in figure 2. Good capture characteristics can be obtained over the Mach number 
range with the fixed geometry design; and spillage, produced by the flow being turned 
toward the opening in front of the cowl by a transient shock system, permits startiig at 
a low Mach number. The downflow produced by the sweep during normal operation tends 
to reduce the static -pressure rise near the top surface and should make possible the inges- 
tion of the vehicle -forebody boundary layer without separation. 

Inlet Startirg 

Inlet starting at the low end of the Mach number range is primarily a function of 
contraction ratio, which is influenced by the amount of sweep, the strut design, and the 
cowl leading -edge location. From simple one -dimensional considerations at an entrance 
Mach number of 3.0, an area contraction of less than 30 percent is necessary downstream 
of the plane of closure corresponding to the cowl leading edge. However, in figure 3 the 
normal plane B-B havirg the minimum cross-sectional flow area is shown to be located 
at the cowl leading edge; therefore, on the average there is no contraction downstream of 
the cowl leading edge. This is a strong indication that the inlet will have no starting prob- 
lem, and early investigations of similar desgns substantiate this conclusion. 

The frontal height-width ratio also contributes to the starting characteristics of 
the inlet. If the width is greater than the height, the inlet is longer and end effects from 
the top surface and cowl begin to dominate the throat. If the width is much less than 
the heght, the struts become slender and structural problems can appear. Based on 
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preliminary investigations, a width-height ratio of 0.8 was selected for the inlet design 
to reduce end effects and to permit a reasonable flight weight structure. 

INLET MODEL DESIGN 
Sweep Angle 

In order to develop the swept shock-wave system for these inlet configurations, it 
was necessary to understand the process that the supersonic flow undergoes as it strikes 
a swept wedge. A detailed discussion of the calculation procedure for obtaining the shock- 
wave ai^le Cjjz) downflow 6jjy, and the flow properties behind the shock wave is given in 
reference 9, and a method for predicting the complete shock-wave train for the inlet is 
developed in appendix A. 

A shock wave may be attached or detached dependii^ on the sweep angle and Mach 
number; consequently, the sweep angle determines the lowest Mach number at which the 
shock waves can be attached at the strut leading edges. If the sweep is too high, low Mach 
number operation will result in shock waves being detached well upstream in the inlet 
ahead of the struts. These detached shock waves may create a situation where disturb- 
ances originating in the combustor may extend upstream of the struts locally and produce 
an undesirable inlet combustor interaction. Too little sweep means the internal contrac- 
tion is high and no mechanism is provided for sufficient flow spillage for inlet startii^ at 
low Mach numbers. Because of the difficulty in analyzing inlet performance at low Mach 
numbers when the shock waves become detached, two models with struts were built and 
tested over the Mach number range of 2.3 to 6.0. These models had sweep angles of 60° 
and 56°. Several smaller models without struts but with sweeps of 50° and 0° were also 
tested at Mach 4. The results of these tests indicated that a sweep ai^le lower than 50° 
could be obtained with adequate starting capability and a sweep angle of 48° was selected. 

Mach 6 Shock-Wave System 

The theoretical shock diagram for a Mach number of 6 is shown in figure 4 with 
tables for the various flow passages and struts within the inlet. Section A-A is a hori- 
zontal plane parallel to the vehicle underbody and cowl. The sidewall compression angle 
was kept low to prevent the possibility of boundary -layer separation due to boundary- 
layer— shock interactions. (See ref. 10.) A detail of the predicted shock-wave structure 
in the vicinity of the struts is given in figure 5, and the properties of each numbered bay 
are listed in table I. The struts provide approximately 75 per cent of the inlet static - 
pressure rise (46-percent decrease in throat area measured in the xz-plane). The side 
struts were positioned within the inlet so that two-thirds of the flow area available was 
in front of the center passages. This flow division permits the fuel to be Injected equally 
from the surfaces of both the center strut and side struts (only six injection planes). 


8 



Boundary -layer transition on the sidewalls is expected to occur ahead of the struts, 
and the correction of 0.4° (based on flat-plate predictions) in conjunction with the side- 
wall angle provides a nominal turning of 6°. Because the chords of the struts are small, 
the thin boundary layer of the struts is expected to be largely laminar in front of the 
throats, and separation of this thin boundary layer is expected to be of little consequence. 
Although it is impossible to prevent the sidewall shocks from merging with the strut 
shocks at some flight Mach number, the situation is relieved for the Mach 6 shock-wave 
system by changing the wall slope of the struts at appropriate locations and either can- 
celing or reducii^ the shock strength. 

Inlet Design Performance 

As previously mentioned, cooling requirements are a major consideration in com- 
bustor design and tend to limit the inlet contraction ratio. The inviscid design aerody- 
namic contraction ratio (which is computed from Mach number at the face of the inlet, 
throat Mach number, and total-pressure recovery) for the inlet at Mach 6 was 6.85 for 
the side passage (bay 6) and 7.57 for the center passage (bay 10) and yielded a mass flow 
weighted average theoretical contraction of 7.3. The geometric throat gaps between the 
struts are greater than those for a two-dimensional inlet of the same contraction, because 
a portion of the contraction is produced when the flow is turned toward the cowl (^xy in 
table l). In this inlet design the ratio of throat gap to inlet height is 0.042 for each cen- 
ter passage and 0.041 for each boundary -layer -corrected side passage. The width ratio 
measured in the xz-plane sidewall leading edges to the throats is 5.74 for the inviscid 
side passages and 6.17 for the center passages when the predicted flow split of 67 per- 
cent for the center passage and 33 percent for the side passage is assumed. 

In the inviscid flow the average throat Mach number is 3.4 and the total-pressure 
recovery is 0.88; as a result, there is an adiabatic kinetic energy efficiency of 99.5 per- 
cent. The inclusion of the estimated boundary-layer losses on the struts, sidewalls, and 
top surface reduces the total-pressure recovery to approximately 0,67 = 98.3 percent^; 

but this value still does not include corner effects. Inlet capture at Mach 6 was predicted 
to be 93 percent for the cowl location discussed in the next section. This predicted cap- 
ture is based on the spillage generated when the flow is turned toward the opening in front 
of the cowl which is computed from matching pressure and flow direction between the 
external and internal streams. 


End Effects 

Although no attempt was made to analyze viscous corner boundary -layer interaction 
regions in the inlet throat, inviscid calculations of the flow near the top surface and cowl 
were made. Because the flow is turned away from the top surface as it proceeds through 
the inlet, a fillet was added as shown in figure 6(a). To match precisely the downflow, the 


9 



top-surfaee contour must vary with local Mach number, but to avoid this complication, a 
single contour angle and position was selected. The angle was kept small {4P) to avoid 
shock interactions which could lead to separation of the thick, top-surface boundary layer. 
The leading edge of the top- surface fillet coincides with the location of the sidewall shocks 
with a Mach number of 5 in front of the inlet. With Mach 6 in front of the inlet, the top- 
surface fillet is upstream of the sidewall shock wave, and a 4° shock wave is produced as 
illustrated in section B-B of figure 6(a). At the throat the fillet displaces approximately 
the same amount or cross-sectional area as it would have if it had been on design for each 
bay, the Mach number in front of the inlet being equal to 6. The corrected values of flow 
parameters near the top surface for this constant-angle fillet are given in table n(a), where 
the angle of flow in each bay has been corrected to match the 4° slope. 

When the internal flow strikes the cowl, the flow must be turned back parallel to the 
cowl internal surface; as a result, a cowl shock wave and a high-pressure region are pro- 
duced. The results of this flow turning are given in table 11(b) where 6^ = 0°. To com- 
pensate partially for this high pressure, the throat area next to the cowl was opened by 
relieving the struts and sidewalls. (See fig. 6(b).) This relief area begins where the 
shock wave from the cowl leading edge strikes the struts and sidewalls, and because the 
center strut is located downstream of the cowl leading edge, the strut does not extend to 
the cowl surface. For the inlet model, this cantilevered center strut was secured to the 
cowl with a pin located behind the strut maximum thickness. 

Selection of the proper cowl leadii^-edge shape was based upon previous model 
testing. Cowls which enclosed the area behind the sidewall shock waves providing 
100-percent capture made the model difficult to start, produced a cowl shock covering 
most of the inlet throat, and generated a very large corner interaction region between 
the sidewalls and cowl. A partial cowl whose leading edges were swept back at 50° 

(f^. 6(c)) was selected, which provided an open area near the inlet throats. This open 
area permits spillage for low Mach number inlet startii^ and provides for some sidewall 
boundary-layer bleed. The exact location of the cowl relative to the struts was deter- 
mined from additional model testii^ as described in a later section. 

Off-Design Performance 

Although the inlet is not a point design at a local Mach 6 in the conventional sense, 
it is designed to give the highest relative performance at this Mach number. As the Mach 
number is decreased as shown in figure 7, the shock-wave system shifts forward, and 
experimental observations on previous models indicate an increased spillage and lower 
aerodynamic contraction. For this design at a Mach number of about 3.5, the shock waves 
become detached from the strut leading edges. With this swept inlet concept, the aero- 
dynamic contraction will increase with increasing Mach number until the shock waves are 
fully attached in the throat region. At still higher Mach numbers the number of shock 
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waves in the inlet decreases as they pass through the throat and the contraction begins to 
decrease slightly. 

The combination of detached shock waves, with resultii^ spillage and variable con- 
traction with Mach number, permits the inlet to have fixed geometry, to start at a low 
Mach number, and to provide enough contraction for successful operation at a high Mach 
number. 


APPARATUS AND TEST PROCEDURE 
The Inlet Model 

A photograph of the inlet model designed with regard to the aforementioned concepts 
appears in figure 8. One sidewall has been removed in figure 8(a) and several pressure 
rakes can be seen. The model is shown upside down, and the 45.72-cm (18-in.) plate 
extendii^ ahead of the sidewalls generates a simulated vehicle -forebody boundary layer. 
The boundary layer from this plate encounters early transition by trips located near the 
leading edge and the resulting boundary -layer profile entering the inlet is measured by a 
three-prong adjustable rake. 

The model is 90.2 cm (35.5 in.) long not including the foreplate, and inlet frontal 
dimensions are 19.05 cm (7.5 in.) high by 15,24 cm (6.0 in.) wide. The aluminum top 
surface was machined in one piece, the forepiate being detachable. Assembly consisted 
of pinning and boltii^ the aluminum sidewalls to the top surface and then the stainless^ 
steel cowl to the sidewalls in any one of three possible positions. Three stainless -steel 
struts were bolted in slots machined in the top surface which were then sealed. The 
struts were attached to the cowl by pins instead of bolts to reduce thermal stresses cre- 
ated by changes in the lei^ths of the struts and heights of the sidewatls. The normal to 
the leading-edge radii of the sidewalls, cowl, and struts was about 0.01 cm (0.004 in.), 
whereas the foreplate leading-edge radius was 0.06 cm (0.023 in.). Partially visible in 
the upper left of figure 8(b) is the actuator mechanism which moved pressure survey rakes 
inside the model. 

Schematic drawings of the model are given in figure 9, Stainless -steel cheeks 
attached to the exterior of the sidewalls simulated adjacent inlet modules up to the inlet 
closeoff station next to the cowl. A survey station for measuring inlet capture was pro- 
vided by a swept flat section of sidewall located downstream of the struts, A sidewall 
relief area next to the cowl is also illustrated, along with the footprint of the side struts 
on the cowl in section D-D. Detailed strut dimensions (fig. 10) and the relative positions 
of the struts and cowl (fig. 11) are measured in any xz -plane parallel to the foreplate and 
away from the relieved area near the cowl. The "station" positions are relative to the 
sidewall leading edges in the same xz-plane. A second center strut, which provides an 
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increase in contraction, is also described in f%ure 10. The difference in contour for this 
second center strut is similarly shown by the dashed lines on some subsequent figures in 
the report. 


Boimdary -Layer Trips 

To insure a turbulent boundary layer entering the throat passages, boundary -layer 
trips were attached to the sidewalls as well as to the foreplate. (See fig. 12.) The diam- 
eter of the steel balls (0.159 cm) was approximately equal to the estimated, flat -plate, 
boundary -layer thickness measured 7.6 cm downstream from the foreplate leading edge. 
The balls were spotwelded to steel strips which were in turn fastened with epoxy to the 
aluminum surfaces. The trips on the foreplate were utilized to increase the thickness 
of the boundary layer which would be entering the inlet, to determine whether there 
were any adverse effects associated with ingesting the boundary layer of a vehicle 
forebody. 


Model Instrumentation 

Figure 13 locates the 116 static -pressure orifices distributed throughout the inlet 
model. The orifices were strategically located to determine inlet starting, pinpoint shock- 
wave position, and aid in evaluating inlet contraction. Because of the variety of locations, 
the position reference varies for each group of orifices. Iron-constantan thermocouples 
were installed in the right sidewall as shown in figure 14, Each thermocouple lead was 
spotwelded to the aluminum surface instead of the leads being welded together, in order 
to determine more precisely the surface temperature. 

Pressure Survey Rakes 

The three-pronged foreplate boundary -layer probe is described in f%ure 15. This 
probe, alined with the sidewall leading edge at Y/H - 0 and Z/H = 0.133, was adjusted 
in height between tests to obtain detailed inlet entrance conditions near the top surface. 

The remaining survey probes were positioned laterally, by an electric motor and actuator 
attached to the model (fig. 16), in one of five access locations. Locations 1 to 4 provided 
for probe surveys across the inlet's throats, whereas location 5 provided access to the 
capture measurement station downstream of the struts. Tubing for the throat survey 
probes (fig. 17) was routed through the hollow actuator shaft, whereas capture measure- 
ment probe (fig. 18) tubes were carried out the rear of the model as illustrated in fig- 
ure 8(a). The throat pitot probes were designed to survey one side and one center pas- 
sage simultaneously, and the static survey probe surveyed only one center passage. It 
was necessary to rely on wall static data for the side passage. The captxire measurement 
station 6.35 cm (2.5 in.) downstream from the struts was surveyed by the seven-prong 
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pitot and static probes of f^ure 18(a). A single, stationary tube (0.102 cm LD.) 
extended through the top surface and was bent toward the flow to obtain pitot pressure 
data in the front of a center passage near the top surface. This tube was bent to differ- 
ent locations across the passage between test runs. 

Because of the static -pressure gradients in the small throat area to be surveyed, 
a conventional static-pressure probe, with the orifices located 10 to 20 diameters down- 
stream from the tip, was found to be xmsatisfactory. Therefore, a new static probe 
design (ref. 11), with the orifices approximately 3 tube diameters downstream from the 
probe tip and on a 3° conical shoulder, was used for both the throat and capture measure- 
ment surveys. These static-pressure probes were calibrated at a Mach number of 4.0, 
and the recorded pressures were found to be in error by less than 5 percent. 

Mach 6 Test Facility 

Figure 19 is a sketch of the Langley 20-inch Mach 6 tunnel. Tunnel test-section char- 
acteristics and flow calibrations can be found in the appendix of reference 12. The tunnel 
total temperature and pressure were normally 467 K and 11.9 atm, respectively, and pro- 
vided a Reynolds number per meter of 9.8 X 10®. One run was made at a reduced Reynolds 
number per meter of 3.3 X 10® at a pressure of 4.4 atm, and all runs in the blowdown tun- 
nel were restricted to less than 2 min. The model was mounted upside down in the center 
of the 50.8-cm-square test section with two 15.24-cm steel channels bolted to the tunnel 
floor as shown in figure 20. 


Facility Instrumentation 

Tunnel pressure was recorded with strain-gage pressure transducers, and model 
static pressures were divided between six 48 -port scanivalves. Pitot pressures and all 
survey data were measured by either strain-gage pressure transducers or multirai^e 
capacitance -type pressure transducers. Pitot position was determined with an electronic 
bridge circuit, and all data were processed by an electronic data processing system.. 



valve stepping mechanism was used to trigger the recording system, once each second, 
while the survey probe was moved continuously across the flow. Before the test the 
speed of the throat probes was selected by varying the voltage to the dc motor until the 
probe would span the flow in one test run. At the capture measurement station two test 
rims were required to span the flow area with the probe. An analysis of the survey data 
also indicated that pressure lag was not significant even thorgh the connecting tubing was 
up to “3 m in length. 
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DATA-REDUCTION PROCEDURE 


A curve-fitting interpolation procedure was utilized to expand the pitot and static 
survey data into a grid network. Mach number, total pressure, and unit mass flow were 
calculated for each grid point; and contour maps of each parameter were plotted by the 
computer's graphic system. Inserted into the program for each grid point was an upper 
limit on total-pressure recovery which was obtained from the inviscid shock diagram. 

(See fig. 5.) If the total pressure was greater than this limit, the measured pitot pressure 
and the limiting total pressure were used to compute the flow parameters including the 
static pressure. The measured static pressure was discarded when the recovery limit 
was exceeded for any particular grid point because of the relative inaccuracy of the static- 
pressure measurements. After completing the grid, numerical integrations were per- 
formed to compute a mass-weighted Mach number and total-pressure recovery for the 
inlet throats and a value for a capture parameter at the capture measurement 

station. 


RESULTS AND DISCUSSION 


Initial Inlet Tests 

The results of the first inlet tests at Mach 6 indicated the model was not operating 
as expected. Initial testing of the model at Mach 6 indicated too much compression and 
possible choking was occurring within the center passage as indicated in figure 21. Upon 
investigation, it was discovered that the 6° sidewall shock wave was striking the side strut 
too near the leading edge, This shock wave then combined with the 4° side-strut shock 
wave and produced a 10° wave which reflected between the side and center struts (fig. 22) 
unlike the expected pattern of figure 5. The location of this sidewall shock wave was also 
determined by removing the struts and observir^ the sidewall static -pressure distribution 
(fig. 23). The error in shock-wave position as measured in the xz-plane was only about 
0.64 cm (A/H a 0.033) measured approximately 36 cm from the sidewall leading edge. 

The effect of the increased compression of the center passage extended across the 
inlet in the vicinity of the cowl to the sidewall, as indicated by the disturbance in the oil- 
streak photograph of figure 24(a). More detail concerning the oil study is provided in a 
later section, but proof that the sidewall disturbance was created by the center passage is 
illustrated in figure 24(b) with the center strut removed. This photograph shows no dis- 
turbance next to the cowl. The measured captvire was only 81 percent with the choked 
cowl compared with 92 percent >Vith the center strut removed. 
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Modified Inlet Conf^viration 


To deal with the miscalculated sidewall shock-wave location, the strut arrangement 
was altered to move the shock waves toward the center -passage throat. The three struts 
and cowl were moved forward AX/H = 0.213 (AX =» 4.06 cm); and the side struts were 
moved toward the sidewalls AZ/lEl = 0.022 (AZ = 0.43 cm) to maintain the same contrac- 
tion and percent of flow in the side passage. The resulting configuration is shown in fig- 
ure 25. It was also observed that the experimentally determined location of the sidewall 
shock waves could be duplicated theoretically by the addition of 0.83° to the sidewall com- 
pression angle. This correction is necessary because of end effects from the top surface, 
or model misalinement, or the use of flat-plate boundary -layer calculations on the swept 
sidewalls and struts. There was some concern for the inlet operation at lower Mach num- 
bers when the shock wave moved forward along the side strut. However, the shock waves 
become detached from the swept compression surfaces and should prevent boundary -layer 
separation and choking by spreading the static-pressure rise along the strut surfaces. In 
fact, unpublished data from low Mach number tests support this conclusion. 

The shock waves were recomputed with the new strut locations and the corrected 
inviscid sidewall compression angle, and the results are given in figures 26, 27, and in 
tables in and IV. Compression in the side passage increased because of the stronger 
sidewall shock wave. Mach number, recovery, and aerodynamic contraction changed 
from 3.51, 87.9 percent, and 6.85 to 3.45, 86.2 percent, and 7.07, respectively. The 
center -passage throat is now divided between bays 10, 14, and 16. The predicted, invis- 
cid, mass-weighted average Mach number, recovery, and aerodynamic contraction for 
the inlet changed from 3.44, 88 percent, and 7.3 to 3.37, 85 percent, and 7.6, respectively, 
for the modified design. 

From the initial inlet tests, moving the cowl to the forward position (fig. 11) aggra- 
vated the choking situation next to the cowl; and the cowl at the most rearward position 
failed to help the design. Therefore, the center position of the cowl relative to the struts 
(fig. 25) was maintained for the new configuration. This movement of the cowl forward 
(AX/H = 0.213) relative to the sidewalls increased the theoretical capture by 2 percent 
to 95 percent. 

The center-passage-throat gap was increased when the struts were moved forward, 
ilthough the inlet aerodynamic contraction increased slightly. The purpose of center 
strut 2 was to decrease again the center -passage -throat gap and to determine operating 
sensitivity of the inlets on this parameter. 
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Inlet Entrance Conditions 

The foreplate static-pressure distribution and the Mach number profile for the flow 
entering the inlet are given in figure 28. The static pressure on the foreplate was above 
the free- stream value because no boundary -layer correction was applied; in addition, the 
boundary layer was thicker than predicted by the flat-plate calculations probably as a 
result of trip losses. 


Wall Static -Pressure Distribution 

Figures 29 to 41 present the results of the static -pressure data throughout the 
model, which are compared with the predicted results of figures 26 and 27, and of 
tables III and IV. The solid symbols are data from orifices used to check flow symme- 
try with the two passages on the right side of the model. The round symbols are data 
obtained with the initial configuration and the x-position of this data has been shifted to 
correspond to the new strut locations. These data indicate in figures 29 to 33 that there 
was reasonable agreement between predicted and measured operation of the side pas- 
sage away from the cowl. Some additional compression was observed on the cowl of the 
modified configuration (fig, 32) that was probably due to corner interaction phenomenon 
and blunt leading-edge effects. The static -pressure distribution along the side -passage 
throat is summarized in figure 33. The symbols of figure 33(b) were obtained from ori- 
fices located on the side strut; and, because no static surveys were made in the side 
passage, a linear pressure distribution was assumed across the throat between the strut 
and sidewall for data analysis. The data from the sidewall throat (squares in fig, 33(a)) 
of the new configuration were neglected because these orifices were no longer at the 
throat but downstream of the struts. 

Data in the center passage are given in figures 34 to 41. The center passage of the 
initial design also operated as anticipated near the top surface (fig. 34), unlike the evi- 
dence of too much compression on the side strut at Y/H - 0.43 (fig. 21 or 35). The 
cowl pressures were greatly reduced with the modified configuration (figs. 38 and 39), 
but no theoretical value is shown because of the complicated flow generated by the strut 
relief (fig. 6(b)). The hump in the center -passage -throat pressure distribution (fi^. 40) 
is attributed to a corner effect from the top surface creating a relocation of the strut 
shock waves. In general the center -passage compression was slightly greater than that 
predicted because no boundary -layer correction was made on the struts. The diamond 
symbols of figures 33 to 41 correspond to the larger center strut. From figure 27 it is 
clear that the larger center strut reduces the gap between the struts and moves the throat 
downstream. None of the experimental measurements would be expected to be affected 
by the larger strut except wall statics on the side strut in the downstream portion of 
bay 16 (fig. 27). Figure 35 does, in fact, show a higher pressure on the most downstream 


16 



static orifice. In all other instances the larger strut had a negligible effect on the 
pressure measurements and it is concluded that the strut could be used successfully at 
Mach 6 to provide added contraction ratio. The static pressure on the small center 
strut (fig. 41) is lower because the orifices were located downstream of the strut shoul- 
der. Appendix B presents pressure levels found within the model for a low-pressure 
test (R * 3.3 X 10® per meter) and also an unstarted condition. 

Oil-Flow Study 

Although static pressures could be, and were, monitored to detect inlet starting, the 
most rapid and reliable method was to observe, by closed circuit television, the oil-flow 
pattern formed as the blackened oil droplets moved across the model's surface under the 
approachii^ flow. If the model did not start, all the oil on the sidewall moved in a 
curved path toward the bottom of the model. The only evidence of any shock wave was 
observed well upstream of the sidewall leading edges on the top surface. With the started 
inlet (fig. 42), the oil on the sidewall moves toward the throat along a line parallel to the 
line of intersection between the top surface and the sidewall. The exception was the oil 
near the bottom edge of the sidewall which followed the path of the flow spilled from the 
inlet. When the oil reached the h^h -pressure region formed by the side strut shock wave 
on the sidewall, it turned toward the cowl along a line approximately parallel to the swept 
leading edges. This oil accumulation line extended out the bottom of the inlet in front of 
the cowl with no interruptions, unlike the original strut configuration (fig. 24(a)). 

Wall Surface Temperature and Boundary - Layer Analysis 

The inlet did not reach an equilibrium temperature (fig. 43) because of a test time 
limit; however, the temperature was fairly uniform because of the high thermal conduc- 
tivity of the aluminum. The dashed lines of figure 43 represent the wall temperature 
distribution selected for the boundary -layer analysis conducted with a modified version 
of the boundary -layer program of reference 13, This integral method computer program 
provided a viscous correction to the top-surface sidewalls and struts and reduced the 
total-pressure recovery from 86.2 percent to 54 percent for the side passage. The center- 
passage recovery changed from 84 percent to 73 percent, and the inlet average recovery 
was decreased from 85 percent to 67 percent. Because the surface of the cowl was small, 
it was neglected in the analysis, but the computed boundary -layer thickness for the rest 
of the inlet was nondimensionalized by the inlet height and plotted in figure 44. Transi- 
tion was forced on the top surface 7.62 cm from the foreplate leading edge because 
of the trips. Natural transition was assumed at s'/H ~ 0.8 for the sidewall where the 
Reynolds number based on momentum thickness was equal to 1000. The sidewall boundary - 
layer trips had been removed prior to the modified configuration test when it was observed 
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that they had little effect on the side-i»ssage experimental data. The boundary-layer 
flow over the struts was assumed to be turbulent from the swept leading edges. 


Throat Surveys 

Pitot pressure distributions for both the side and center passage are presented in 
f%ure 45, where W is the throat gap and Z' is equal to zero at the sidewall for the 
side -passage throat and equal to zero at the side strut for the center -passage throat. At 
the top of the model (Y/H = 0.14), the center throat data were obtained from the fixed 
tube, which was bent to a new Z' position for each test. The side-passage pitot profile 
(Y/H = 0.17) was obtained from the theoretical boundary -layer profile at Z/W - 0.5 
and then assumed similarity with the profile at the adjacent side throat station of 
Y/H = 0.26. The solid symbols were not data points but depended on theoretical boundary 
layer calculations to extend each survey (dashed lines) to the wall static value for data 
reduction. 

Data station Y/H = 0.43 was considered the station least influenced by top-surface 
and cowl flow effects, and surveys for the low Reynolds number test and the large center 
strut test are presented for this station in f%ures 45(g) and 45(h). The dashed lines in 
these two figures were taken from the standard tests (fig. 45(c)) and indicate that neither 
the large center strut nor the low Reynolds number tests had any significant effect on 
inlet operation. 

The static survey data (fig. 46), obtained for the center passage, were faired to the 
wall values as indicated by the solid symbols. As mentioned earlier, when the measured 
static pressure was low and the total-pressure recovery exceeded the set upper limit 
(88 percent), a new static pressure based on the upper recovery limit was computed as 
illustrated by the dashed lines. To check the effect of this restriction, it was determined 
that a 10-percent increase in the limit increased the mass-weighted total-pressure recoV' 
ery by approximately 3 percent. The straight lines are the static -pressure distributions 
assumed for the side passage and are based on ne^hboring wall static values. Other 
static -pressure data estimated from ne^hborhood statics for both the center- and side- 
passage throats near the top surface and cowl are given in figure 47. These additional 
estimated distributions made possible an analysis of a greater jmrt of the throat flow ares 
as defined by the pitot surveys. 

Jhe throat stations at Y/H = 0.43 were selected to be compared with the theoreti- 
cal inviscid and boundary -layer calculations, and the resultii^ Mach number profiles are 
found in figure 48. The agreement is good for both passages; the dip next to the center- 
passage side strut naay be caused by a probe-tip shock-wave interaction. The small hash 
marks indicate the pitot measurement closest to the sidewall or strut, and the remainder 
of the curve is the result of fairing the pitot pressure to the wall static -pressure level. 
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Throat Contour Maps 


Contour maps of the results of compiling the throat pitot and static data are given in 
figures 49 to 52. Figures 49 and 50 are maps of the data input used to compute the Mach 
number and recovery maps of figures 51 and 52. Each map is shown with the width scale 
seven times the height, which makes the relief area next to the cowl (illustrated in fig. 51) 
appear to be out of proportion. The side-passage Mach number contours (fig. 51) are rela- 
tively symmetrical, and the prediction of boundary -layer thickness 6 agrees well for the 
top and side surface. The top-surface boundary -layer thickness for the side passage had 
to agree because of the imposed boundary conditions. A nearly horizontal shock of about 8° 
turnii^ was generated by the cowl leading edge and is still near the cowl surface at the 
throat. This discrete shock is smeared by the computer programs interpolation process; 
and consequently a vertical Mach number gradient extending well beyond the predicted 6 
for the cowl is indicated. There is some roundii^ of the contours at the corners, but no 
flow separation is detected. The Mach number contours for the center passage are not 
as symmetrical because of the greater shock-wave concentration (fig. 51(b)); however, 
the Mach 3 (mass-weighted average equals 3,11) contour encloses the major portion of 
the total area. The mass-weighted total-pressure recovery recorded in figure 52 for 
both the throat passages was obtained by a computer program which averaged the values 
of approximately 1000 grid points, evenly spaced over the throat area. One case was 
also integrated by hand with negligible difference, and the results verified the computed 
results. The mass-weighted average recovery for the two passages is 59 percent when 
the losses on the foreplate are included and 61 percent when they are neglected. The 
central area of nearly constant total-pressure recovery in f%ure 52(b) is in part a result 
of the assumption of a total-pressure recovery limit in regions where the measured 
static-pressure level was too low (fig. 46); however, this assumption is considered to be 
justified on a phenomenological basis as well as by the similarity in shape between the 
measured and derived static-pressure profiles (fig. 46) and the reasonable agreement 
between theory and data (for example, fig. 48). 


Capture Measurement Results 

The procedure for analyzing the flow at the capture measurement station was iden- 
tical to that of the inlet throats. At this station static survey data were taken at each pitot 
survey location. Each rake had seven probe tips. The data from which the Mach number 
and capture parameter pv/pjVj were derived for figures 53 and 54 are discussed in 
appendix C. The wakes of the three struts, which are about 7.5 cm upstream, are 
detectable in the Mach number 'map (fig. 53); in general, lines of constant Mach num- 
ber are parallel to the sidewalls. Besides measuring inlet capture flow, the capture 
parameter (fig. 54) is a good indicator of flow gradient direction because it is less sen- 
sitive to static -pressure error than either Mach number or recovery. The average 
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value of P^/Pi^i was 3.18 for the inlet at the capture measurement station; before 
capture flow could be calculated, however, some estimate of flow direction to de^te^^ 
cross-sectional area had to be made. The assumptions in figure 55 are: flow parallel 
to the top surface in the top -surface boundary layer; flow parallel to the cowl below the 
estimated location of the cowl shock wave; and flow down at 8° for the remainder of the 
area. This 8° downflow was the flow turnii^ which was computed from the rise in static 
pressure due to the cowl shock. The flow was also assumed to be parallel to the side- 
walls at this station. With these restrictions a capture of about 94 percent was computed 
for the inlet at Mach 6. 


Performance Results 

Because the struts were positioned within the inlet to provide two-thirds of the flow 
to the center passage, the captured flow (94 percent) was assumed to be split between the 
center and side passages in the ratio of 63/31. With this criteria the Mach 6. integrated 
performance parameters based on the tunnel free -stream conditions are tabulated for 
each passage, and the total inlet, in figure 56. The side-passage total-pressure recovery 
was lower than the center -passage recovery because of the relatively thicker boundary 
layer on the sidewalls. The average viscous total -pressure recovery was 0.59 compared 
with the predicted value of 0.67 for the initial inlet conf juration, which did not include 
corner effects. The aerodynamic contraction ratio, which is based on the average throat 
Mach number and total pressure, is 7.0 instead of the design value of 7.3. This is because 
moving the struts upstream increased the throat width somewhat and because the meas- 
ured total-pressure recovery was slightly lower. The larger center strut increased the 
contraction, but no data were taken at the new center -passage throat location. Included in 
figure 56 are curves from reference 7 predicting the inlet kinetic energy efficiency and 
capture over the flight Mach number range of 4 to 10. The Mach 6.0 inlet data have been 
entered in the figtire for a flight Mach number of 7.6, which indicates a representative 
amount of vehicle forebody compression. The measured kinetic energy efficiency as 
determined by adiabatic process was 97.7 percent compared with the predicted tunnel 
value of 98.3 percent. The predicted curve is slightly high primarily because the cowl 
shock and viscous corner interactions were not included. The measured captured mass 
flow from in front of the inlet was 94 percent and matched the predicted value. The theory 
is expected to be less accurate at low Mach numbers because of the formation of detached 
shock waves. (See fig. 7.) 


. CONCLUDING REMARKS 

As part of a Lai^ley research and technology program focused on the development 
of a concept for an airframe -integrated scramjet engine (Langley Scram jet Module), a 
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detailed performance evaluation of the baseline inlet configuration at Mach 6 (simulated 
flight Mach number of approximately 7.6) has been conducted. 

Mach number profiles in the inlet throat agreed reasonably well with the predicted 
results both for the inviscid flow and the boundary -layer calculations on the struts and 
sidewall. The mass- weighted average throat Mach number was 3.0 for the side passage 
and 3.1 for the center passage. This value compares with the inviscid values of 3.4 for 
each of the two passages. This additional compression was produced by the boundary 
layer, viscous corner interactions, and other end effects (for example, the internal cowl 
shock). 

An adiabatic kinetic energy efficiency of 97.7 percent (0.59 recovery) was measured 
and compared with a predicted value of 98,3 percent (0.67), which does not include corner 
or end effects. 

The average inlet aerodynamic contraction ratio was 7.0 instead of the predicted 
value of 7.3 because of the slightly lower total-pressure recovery. However, the results 
indicated that the contraction ratio can be increased by the use of a larger center strut, 
which operated successfully but no survey data were taken. 

The measured inlet capture flow was 94 percent which agrees with the predicted 
value of 95 percent. 

The difficulty and importance of properly positioning the shock waves in the throats 
of hypersonic inlets were demonstrated, but no adverse effects were noted as a result of 
the inlet ingesting a boundary layer on the top surface which simulated the boundary layer 
that would be ingested from the vehicle forebody. Overall, the inlet performance is well 
within the acceptable range for high engine performance. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 22, 1975. 
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APPENDIX A 


ANALYTICAL CALCULATIONS OF SWEPT SHOCK WAVES 

Shock-wave systems for swept inlet configurations require a three-dimensional coc 
dinate system to locate the shock waves correctly and to compute the flow properties. A 
selected two-dimensional coordinate system (the xz-plane of fig. 4) is also helpful in mai 
taining visual contact with the problem. Figure 57 illustrates the flow striking a swept 
wedge (surface AGFED) and helps to describe the development of the swept shock wave. 
Points A, H, I, and D lie in the xy -plane and points A, H, G, and B lie in the xz-plane. 

The flow strikes the leading edge at point A; and if the wedge were not swept, the flow 
would only be turned away from the xy-plane by angle and would follow the path AC 

However, sweeping the wedge requires the flow to traverse the surface along AF, the flo 
also being turned away from the xz-plane as shown by angle ^xy which is measured in 
the xy-plane. A swept shock wave (plane ABCD) attached to the leading edge is producec 
and is located with angle which is measured in the xz-plane. As long as the wedg< 

is assumed to be of infinite length, no end effects are encountered; if the wedge is assunc 
to extend from the xz-plane, a reduced pressure, nonuniforni flow region will exist in th< 
proximity of the xz-plane. To eliminate this region, a fillet (AGFB) is added which 
extends out to the shock wave and fills the void left by the flow beir^ turned away from 
the xz-plane. 

A detailed discussion of the calculation procedure for obtaining the shock wave 
angle 6jjz, downflow 5xy» aiid the flow properties behind the shock wave is given in ref 
erence 9, which describes the flow velocity being broken into vector components normal 
and tangential to the leading edge as shown in the sketch of figure 57. Although the velo 
ity component tangential to the leading edge V^ remains unchanged = Vt^2)> 
component perpendicular to the leading edge Vn^l is reduced because it is turned by 
when it encounters the wedge. The vectors on the wedge surface, ^xi,2’ "^t.2» 

then combined to obtain the velocity and direction of flow on the wedge surface. One lin 
tation to the procedure occurs when either the sweep angle or the wedge angle is too grc 
and the velocity component Vn 2 becomes subsonic. The shock wave may become 
detached from the wedge leading edge, as illustrated in the sketch, and the downstream 
flow is not uniform. 


The Swept Shock-Wave System 

To determine the shock- wave orientation and flow properties for a train of shock 
waves as illustrated in figure 58, additional steps are required to compute the shock anj 
in the xz-plane. The computation must include the change in sweep angle and normal 
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APPENDIX A - Continued 


turning angle as the flow crosses each successive shock wave and is turned further from 
the xz -plane. 

A computer program has been written to compute an inviscid shock wave system for 
an inlet, when the inlet geometry and the initial Mach number are specified. Inlet geom- 
etry consists of plane wedges defined by sweep angle measured in the xy -plane and flow 
turning angles which are measured in the xz- plane. The shock-wave computational pro- 
cedvire is described for the three shock waves of figure 58 . The second shock wave 
reflects from a plane of symmetry aloi^ line AB, strikes the sidewall (line BB), and is 
again reflected. The flow vectors behind the three shock waves are labeled V3, 
and V4, respectively. Because the flow ar^le with respect to the top surface (xz-plane) 
increases as the flow crosses each successive shock wave, a new fillet is required to 
eliminate three-dimensional end effects. The downstream fillets have complex orienta- 
tions which are functions of Mach number; but, because they are small and considered to 
have minimal influence on the flow, these effects were neglected in the inlet design. 

The sidewall is the generator of the first shock wave and may be considered to be 
a wedge or wing with sweep Aj measured in the xy-plane, angle of attack and 

dihedral as measured by the angle (i> j. The second shock wave is reflected from the 
plane of symmetry along line AA and is illustrated in detail in figure 59 . For this 
reflected wave the computer program solves the problem of flow across a wedge or wing 
with sweep A2, angle of attack zero dihedral (<^2 “ reference axes 

are x’, y’, and z’. Because the flow has been turned toward the y-axis by the first 
shock wave, the sweep angle A2 increased from the value of Aj. Once V3 and the 
shock-wave angle e^xz)’ been computed, the shock wave is defined in the original 

xz-plane with angle 

The third shock wave is illustrated in detail in figure 60 , where the flow in front of 
the wave V3 approaches a swept wedge (sidewall) which is in the x", y", and z” coor- 
dinate system. The sweep angle is A3. The angle of attack is dihe- 

dral is <t>2‘ The computer program treats this wave in the same manner as the first 
wave to compute £(xz)"» then defines the shock wave in the xz-plane with 

The first and third shock waves are coded "type A" waves by the program whereas 
waves reflected from the plane of symmetry are coded "type B." The flow turning across 
each wave, measured in the xz-plane, can be put into the program, which means the 
strength of the reflected wave (type B) may be made different from the strength of the 
incident (type A) wave. To compute a shock train, the program always begins with a wave 
of type A, but the wave types do not have to alternate as the wave type of each shock wave 
is input to the program. When two type B waves are together, the nomenclature must be 
reversed and the second type B wave considered to be type A; then this new orientation is 
continued. 
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APPENDIX A - Continued 


The program transfers all shock angles back to the xz-plane, where the angle is 
measured with respect to the flow direction in front of the wave for a type A wave; and 
the shock angle is measured with respect to the flow direction behind the wave for a 
type B wave. For a selected path of shock waves the sweep angle and dihedral for the 
first wave are input to the program. Subsequent sweep armies are internally computed, 
as is the dihedral, which is a function of the wave type. The number of shock waves in 
the train is input, and the program will continue to calculate across successive shock 
waves until shock-wave detachment occurs. 

Because the flow properties in oblique shock- wave systems are path dependent, it 
may become necessary to iterate on pressure and flow direction, when waves of different 
turning strengths are encoimtered. Usually, however, the differences in shock-wave 
turning angles are small enough to insure that such effects can be neglected, if a correc- 
tion is deemed to be necessary, an iteration can be done either by hand as the program 
computes across one wave at a time, or by followii^ several flow paths with the program 
and averaging the results in the selected downstream flow bay. 

The Computer Program 

The program is written in FORTRAN IV and is adapted to the CDC 6600 computer 
located at the Langley Research Center, Hampton, Virginia, Although the primary pur- 
pose of the program is to compute flow properties for swept, weak, oblique shock-wave 
systems, additional versatility is available as shown in the input listing attached to this 
appendix. Either perfect-gas or thermally perfect-gas (gas with caloric imperfections 
as defined in ref. 14) problems may be computed by the listed program. Also operational, 
but not included in the program listing, is a subroutine for real-air calculations. This 
real-air subroutine uses either thermodynamic tables or equations of air in thermochem- 
ical equilibrum to compute flow properties behind shock waves. 

For the thermally perfect gas, the local total properties are computed; and the 
specific heat ratio y is computed as a function of static temperature from equation (180) 
of reference 14. The gas characteristics that are assumed for air and are used cxirrently 


in the program are: 

Molecular vibrational-energy constant, Sy 3076 K 

Perfect gas specific heat ratio, Vp 1.4 

Gas constant, R. . . . 1545.31 

Molecular weight 28.9644 


Both weak and stroi^ shock waves and perfect-gas Prandtl-Meyer expansion calculations 
are possible. For the expansion fan, the angles of the leading, trailing, and average waves 
are printed. 
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The dihedral (p is measured in a plane perpendicular to the wedge ridge line 
(line AG of fig. 57) and is defined in a parameter ip where 

, _ ^ cot A tan (p 
^ 1 - . - ■- ■ 
sin 6xz 

This parameter is defined in reference 9 as the ratio of the lengths of two line segments 
(DE/AD) and is used to define the location of the wedge leading edge with respect to the 
xy-plane. For the first wave of a shock train, the wedge dihedral ^ does not have to 
be calculated because the leadir^ edge is assumed to lie in the xy-plane and has a 
value of zero. For a single -wave calculation the leading edge may be lifted out of the 
xy-plane; and when the dihedral is zero, ip has a value of 1. 

Program Input 

Card 1 - FORMAT(4F10.4,2F5.2,2F10.4,5I1) 

OC Number of successive waves (Maximum = 20) (cols. 1 to 10) 

XM Initial Mach number, M (cols. II to 20) 

ALP Sweep angle, A, deg (cols. 21 to 30) 

El Wedge ridge angle, 5xz> (cols. 31 to 40) 

DEADR Dihedral parameter, \p = 1 - ^ (cols. 41 to 45) 

sin 6xz 

GAM Specific heat ratio, y, for perfect gas (cols. 46 to 50) 

PI Static pressure, pj, psia (cols. 51 to 60) 

Set equal to 1.0 if left blank 

T1 Static temperature, Tj^, °R (cols, 61 to 70) 

Set equal to 500° R if left blank 

K1 Type of gas calculation (col. 71) 

0 perfect gas 

2 thermally perfect gas 

K2 Flag indicates value of turning ar^le for each shock wave to follow on 

card 2 (col. 72) 

0 all waves will have 6xz * El 

2 OC values of will follow on card 2 

K3 Flag indicates wave type for each shock wave to follow on card 3 (col. 73) 

0 no card 3 necessary, waves will alternate Type A and Type B 
2 OC wave types to follow on card 3 
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K4 Flag indicates strei^th of shock wave (coL 74) 

0 weak shock waves 

2 strong shock wave (limited to OC = 1) 

K5 Special input case for dihedral (col. 75) 

0 DEADR = Dihedral parameter, x{/ 

2 DEADR = Dihedral angle, deg 

Card 2 - FORMAT (16F5.3). Do not use if K2 = 0 

EllA Value of Sjjz wave, deg 

Expansion waves are entered as negative 

Card 3 - FORMAT (16F5.3). Do not use if K3 = 0 

TYPE Wave types 

1. Type A 

-1. Type B 

The first wave of a shock train must be type A. 

A listir^ of the computer program follows: 



PROGRAM inlet ( input , output, tap E5= input,! APE6=0UTPUT , taper ) 

10 0100 

c 



700QQQ 

c 

SWEPT SHOCK WAVE PHUGPAM A19R1 BY C.A. TRFXLER 10/22/74 

300 100 

c 



4QQQQQ 

c 

INPUT FOP 

CARD 1 FnKMAT(4Fn,4,2F5,2,2Fn*^,Em 

500000 

CNO* WAVES MACH 

N0» SWEEP DELTA DIH GAM PI Tl KKKKK 

600000 

C 


DEG, DEG. OE/AO P$IA OEG.R 12345 

700000 

C 



BQOQOQ 

C 

SWEEPS 

SWEEP ANGLE MEASURED IN XY PLANE 

90QOOO...._ 

C 

0ELTA= 

FLOW TURNING ANGLE IN XZ PLANE 

1000000 

C 

DIH = 

dihedral PAKAMETEK (DE/AO) 

1100000 

C 

GAM = 

SPECIFIC HEAT RATIO ** SET EQUAL TO 1.4 IF LEFT BLANK 

1200000 

C 

pi= 

STATIC PRESSURE, PSIA ** SET EQUAL TO 1»0 IF LEFT BLANK 

1300000 

C 

Tl- 

STAT TEMP^, DEG R =<=* SET EQUAL TO 500* IF LEFT BLANK 

1^00000 

C 

AOOITIONAL ENGLISH UNITS 

15COOOO 

C 


DENSITY — LBM/FT**3 

1600000 

C 


VELOCITY - FT/SEC 

1 700000 

c 


ENTHALPY - BTU/L8M 

IBOOOOO 

c 



19C0000 

c 

Kl= BLANK, 

PERFECT GAS /// Kl= 2, THERMALLY PERFECT GAI 

2000000 

c 


Kl= 4, REAL AIR - EQUATIONS (IF Ti = 0 Pl^ ALT I TUOE , . FT* i 

21C0300 

c 


Kl=: 6, REAL AIR ** TABLES (IF T1 = 0 Pl^ ALTITlJDF, FT*L 

2200000 

c 

K2= BLANK, 

WAVE strengths not INPUT /// K2^ 2, WAVE STRENGTHS INPUT 

2300000 

c 

K3= BLANK, 

WAVS TYPES NOT INPUT /// K3= 2, WAVE TYPES INPUT 

2400000 

c 

K4= BLANK, 

WEAK SHOCK WAVE /// K4= 2, STRONG SHOCK WAVE 

2500000 

c 

K5= BLANK, 

OIH = DE/AD /// K5= 2, SPECIAL CASE 

2600000 

c 

^'SPECIAL CASE« DIHEDRAL (OIHJ OEFINEO BY AN ANGLE MEASURED 

2700000 

c 


IN A PLANE PERPENDICULAR TO RIDGE LINE* 

2800000 

c 


SWEEP MEASURED WITH LEADING EDGE IN XY PLANE, 

2 900000 

c 


AND WEDGE AT ZERO ANGLE OF ATTACK. 

3000000 

c 



31QOOOO_ 

c 

INPUT FOR 

CARD 2 F0RMATC16F5.2I NOT INPUT IF K2= BLANK 

3200000 

CDiU 

0(2J DI3) 

0(4)®«*b^D(NW) 

33COOOO 

c 


FLOW TURNING ANGLE IN XZ PLANE (DEG) 

3400000 

c 

NW = 

NO* WAVES 

3500000 

c 



3600000 _ 

c 

INPUT FOR 

CARO 3 F0RMAT(16F5*2) NOT INPUT IF K3= BLANK 

3700000 

CT(l) 

T(2) T(3S 

T(4)*«*«.T(NW) 

3800000 

G 

T= 

WAVE TYPE ««««#* T= 1* FOR TYPE A /// T= -1. FOR TYPE B 

3900000 

C 

NW= 

NO* WAVES 

4000000 

c 



4100000 

000003 

COMMON C, 

THV,GAMP,RBAR«XMWT,Pl,Tl9XM,XMN0RM,T0 
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50003 


COMMOiN P2» T2f OeLTN,GAM,XM2NtEPO,RHrjR 

4300000 

30003 


COMMON EMAX,OELrNM,N,KG*AS,MC 

A 4013 ad 
. . . 4510100 _ 
4600003 _ 
4700000 

D0003 


COMMCN epdF»TPKA» ARfXMlT.XMZT 

30003 


COMMON NOEBUG 

00003 


DIMENSION Ell AL20) . TYPE (20) 

00003 


C=3. 1^159265/180. 3 

4331130 . 

. 49Q31Q1. 

. fllOOlQ . 
3IQa01O_ 
52C33Q0 
5300000 

00005 


GAM^l.400 

00006 

odoor 


TMV=5537,0 

GAMP=1»4 

00010 


K6AR = 15’45.31 

00012 


XMWT=28.Q644 

00013 


NDEBUG=0 

.. 5401000 . 
. 550.0100 _ 

5 600000 

5700000 

0001^ 

00016 

5 

jj=Q .... ... .; . 

K6AD(5f 101 >0C f XM, ALPfEl ,DEA0R,0AM,P1 , T 1 , KGAS , NDF LT , NORD^R , N . KSP 

0005 A 

101 

FORMAT (4Fid.4,2F5.2.2Fl 0.4, 511) 

00054 

00057 

6 

if (EOF, 5 >6, 70 
STOP 

5800000 

5903000 

00061 
00061 
00062 
0006 3 
0036 5 
00071 

jn_ 

CONTINUF 

to = 0. - - - . - - . 

DIH^OEADR ' ’ " 

nstr=n 

IF (GAM.lt. .001 )GAM=l’.4 
IF{KGAS.GT.3)GQ TO 400 

..6111110 

61C0000 

6 ' 2 ododo 

6303310. 

6400000 

00075 

00077 

ooin 

ooioi 

00133 • 

^OQ 

1F(P1.“eq. 3.>P1=U 
IF<Tl.EQ.O.m = 500. 

CONTINUE 

IF(KGAS.Eg.O)KGAS=0 

IF(NUELT.Eg,0)NOELT=0 

6600000 

6700300 

6800000 

6900300 

7oooddo 

00105 


IFtNORDFR .EO.O )N0R6ER=b 

7103000 

00107 


IF (NSTR. eg. 0>NSTR=3 

.720.0000,. 

7300000 

7400000 

00 Hi 
00113 


IF(KSP.e0.0>KSP=0 

IF ( ALP.EQ.O. )G0 TO TOO ' 

00114 


IF(ei.£Q.O.>Gd TO 700 

750 >030 

00115 
QOlW 

700 

IF (KSP.Eg.2#DEADR = l.-TAN{5EADR*C)/(Sr;i(ei«C)<'TAN{ ALP#C> ) 
CONTINUE 

7600000 

7703003 

00136 

00137 


xminf=xm 

PINFK^l.O ■ ■ 

7800000 

■^903000 

ODl^ 1 


TINFR-1,0 

3000000 

00141 


PINFR=1.0 

8133000 

03142 

001^4 


KKGAS=KGAS 

IF(KGAS.Eg.3)KGAS=l 

8200000 

8300030 

^3146 
00150 


MC=bC 

DO 3 I-1,MC 

3_400000 

8503030 

0015 2 


EllAUi^^El 

8600000 

00154 


TYPE(n=0. 

8700000 

00155 

3 

CONTINUE 

8800000 

00157 


IF(NDEIT.EQ.O)GO TO 4 

0900000 _ 

100160 


REA0(5,102H E11A(N0> ,N0 = 1,MC) 

9000000 

00172 


ei=EllA(l> 

9100000 _ 

• 00174 

102 

FORMAT (16F5. 3) 

9200000 

00174 

4 

CONTINUE 

9300000 

100174 


IF<N0RDER.EQ.2iREA0(f ,l02i(TYPE(N0HN0=l.MCI 


00211 


OC=1.0 

9500000 

00212 


1 = 1 

9600000 

00213 


COI3PT=0.0 

9700000 

»00214 


ARCA=1.0 

9813100 

00215 


TPRA=1«0 

9900000 

00216 


K=1 

10000000 

00217 


NN=0 

10100000 

00220 


JJ = JJ<-i 

10200000 

>30222 


AL1=ALP 

10300000 

00224 

7 

CONTINUE 

10400000 

) 00224 


PRINT 106 

10500000 

00230 


FORM'AT(IHI) , 

10600000 

)00230 


PRINT 126, JU,K 

10700000 

00240 

126 

F0RMAT(1X,7HCASE NOi 13 , lOX, 7HW A VE NO, 13, 39X2HK 1, 2X2HK2, 2 X2HK 3, 2X 

10800 300 


-1 

L 2HK4,2X2HK5) 

109QOQQQ 

) 00 240 


1F(NSTR.EQ.2.ANO.E1.GE..01)PRINT 121 ,KKGAS ,NDF LT , NQRDER, NSTR ,K SP 

11000000 

>03267 


IF (NSTRoNE«2*AND«E1*GE*0* JPRINT 122 , KKGAS ,NOELT , NORDEH, NSTR, KSP 

lllOOQQO 

00315 


IF(NSTR.NE.2.AND.E1.LT. J.>PRINT 123 , KKGAS , NOEL T, NORDER , NSTR, KSP 

112C-0000 

00343 


IF(NStR«EQ«2«AN0«El«LT««01 IPRINT 1 24, KKGAS * NDE LT .NORDER. NSTR .K SP 

11 330000 

03372 


lF(KSP.e0.21_PRJNT 125, DIH L 

114.00000 . 

)00402 

121 

F0RMAT(45X14H( STRONG SHOCK* *9X614) 


>00432 

122 

format (45X14H( WEAK SHOCK) ,9X614) 

116Q03Q0 

>00402 

123 

F0RMAT(45X14H( EXPANSION) ,9X614) 

11700000 

>00402 

124 FORMAT (45X14H( NORMAL SHOCK ) .0X6 14 ) 

~ 11 eOOOQQ 

>00402 

JL25 

FaRMAT(70X26H(DIHE0RAL = SPECIAL CASE ^,F7.3,5H DEGl) 

~ I190Q000 

>00432 


_Jri ^ 
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000403 

600 CONTINUE 

„ 1210Q00J)_ 

000403 

ALPH«ALP«C 

12203000_ 

000405 

E=ei*c 


000406 

if(alph,gt®»oougo to 9 

12460^60 

000412 

XMIN-XN 

125000aQ_ 

000412 

XMIT^O. 

12600000 

000413 

XMNORM=XM 

12700000 

000414 

0ELTN=E 

12 80)000 

000416 

OELTNP=^OELTN 

12900000 

000420 

GO TO 14 

13000000. 

000420 

9 CONTINUE 

13100000 

000420 

AC0U*90«*C-ALPH 

13200000^ 

000423 

CD=0C*TAN(AC001 

||[||||||^^ 

000426 

00*0C/C0S(AC0D1 

13400000 

000431 

BC=OC^TAN(EI 

13500000 

000434 

AO=BC 

15600000 

000436 

DE=AO*DEADR 


000437 

PSI=ATAN2IDE,ODi 

13800000 

000442 

OE={OD»*2«-DE**2i**-5 

13900000 

000430 

08=(BC*#2<-0C#*2I*««5 

14000000 

000456 

AE=ABS(AD-DE> 

14100000 

000461 

6E-(CD**2«-AE#*2V#**5 


000467 

THETA=ACOS(OC/OE) 


000473 

ABO£=ACOS( (0B**2<-0E4«2-BE#*2) / ( 2.0*08*0E ) ) 

14400000 

000504 

XM1T=XM*OC/Oe 

i45o_aoao_ 

000507 

XM1N=XM4^S INI THETA * 

14600000. 

000512 

IE(XM1N-1.0J203,203,8 

147G0000 

000515 

8 CONTINUE 

14800000 

000515 

OG=OC*OC/Oe 

14900000 

000517 

CG=UC*SIN( THETA* 

15000000 

000522 

HG=OG=<'T AN( A60E* 


000525 

0H=(HG**2<-0G**2)***5 

15200000 

300532 

HC=(OH**2frOC*+2-2*0*OH*OC#COS(E* 

15300000 

000544 

DeLTN=ACOSUHG4<'24-CG**2-HC**2*/(2,0*HG«CG* J 

15400000 

000554 

IF(eilA(K).LT.O.)DEtTN--OELTN 

15500000 

000560 

deltnp^oeltn 

15600000 

000561 

XMN0RM=XM1N 

15700000 

000563 

14 CONTINUE 

15800000 

000563 

GR=(GAMfl. */iGAM-l*> 

15900000 

000567 

EMAX=SQRT((GAMflV**U.i-(GAM-l. >/2«*XMlN#*2fr(GAM«'l,) /16.*XM1N*<'4) 1 

16000000 

000605 

EMAX=U./(GAM*XM1N**2* **( (GAM4-1® */4«#XMlN^*2«l. fEMAX* 

16100000 

000620 

EMAX=ASIN(SQRT(EMAXI * 

16200000 

000624 

OELTNM =(GAWn« **XMlN=>’«'2/i2«*iXHlN4*2* ( S INC EMAX *«*2 > -!• ) * 

16300000 

000635 

OELTNM= ( DELTNM-1 . * *T AN ( EMAX ) 

16400000 

000642 

DEL TNM^ AT an 11 « /DEL TNM ) 

16500000 

000645 

DLIMIT=DELTNMnO»*C 

16600000 

000650 

IF (XM1N«LT«2» *DLIMIT=^0ELTNM4-2«*C 

16700000 

000656 

IF(DLIMIT*GT«OELTNJGO TO 15 

16800000 

000662 

KGAS=1 


000663 

15 CONTINUE 

17000000 

000663 

lF(KGAS-lJiO,10*ll 

17100000 

000666 

11 CONTINUE 

17200000 

000666 

IF(KGAS«GT«3)CALL RA IR ( ZFT, NSTR* 


000673 


17400000 

000677 


175QQQQQ_ 

000700 

IFCN*1T«100>GO TO 16 


000703 

KGAS=l 

17700000 

000703 

IF(0EtTNH,lT*0ELTN*G0 TO 10 


000706 

IFCN.E0«100*G0 TO 5 

1790QOOO_ 

000707 

16 CONTINUE 

1800 0000 

000707 

TR*T2/T1 

18100000 

000711 

PR1=P2/P1 

18200)00 

000713 

RR=RHQR 

18300000 

000714 

EP^EPO^C 

18400000 

00071 6 

EPF*EP0F4^C 

1850Q.a0Q_ 

000717 

IFCALPH*LT»»001)G0 TO 21 


000722 

10 CONTINUE 

18700)00 

000722 

IFIKGAS-U12»12*25 _ 

18800000 

000725 

12 CAtL PGAS 

18900)00 

000726 

IFIN.EO.IOOJGO TO 5 

19000000 

000730 

TR=T2/T1 * 

19100000 

000732 

PR1=P2/P1 

192QQOQO_ 

000734 

RR=RHQR 


KEf^EiSHIl 

EP=6P0«C 

194QQQX)Q_ 

000737 

EPFt£PDF*C 


000740 

IF(ALPH.GT..001>GO TO 13 

19600000 

000744 

21 CONTINUE 

19 700 OOP 

000744 

ARCA^ARCA*AR 
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000746 

0 ELTN 1 =DELTNP/C 

19900 300 

_000747 

DitTM=X)EtTNH/C 

20000000 


E^IAXl = EMAX/C 

20100300 

000752 

EEF=l.- « 2 , / t GAM- 1 , )*( ( X^TPRAl**( (GAM- 1 . ) /GAM) - 1 . )/XMlNF** 2 ) 


. OOOtAfe 

PiNFR=PINFR=»PRl 

203 GQ 000 

000770 

TINFR=TINFR#TR 

20400000 

000771 

RINFR=RINFR=^RR 

20503030 

000773 

PRINT 115 

20600000 

000 T 77 

115 FORHATl// ^ SUMMAR 7 L 0 F F_ 10 M PROPERTIES ACROSS UNSWEPT WEDGE 

2 Q 703 Q 3 Q 


1 ♦) 

20800000 

000777 

PRINT 17 , XM*ALP*DELTN 1 ,EPD,P 1 ,T 1 

20900000 - 

001017 

17 F 0 RMAT(/lX. 3 HMl=,F 8 . 3 * 7 X 6 HSWEEP=,F 7 - 3 , 7 X 6 HDeLTA=,F 7 . 3 , 12 X 4 HEPS=. 

ZlOOOOdO. 


1 F 7 . 3 . 7 X 3 HP 1 -,F 8 . 4 , 6 X 3 HTI=,F 8 . 2 ) 

21100000 

001017 

PRINT 19 . XM 2 N,PINFR*RINFR.TINFR.GAM__ 

21200000 

001035 

19 F 0 RMAT(/ 1 X, 3 HM 2 =,F 8 . 3 . 5 X. 7 HP/PINF=,F 8 .A. 5 X 7 HR/RINF=.F 8 .A. 3 X 7 HT/TIN 

21300010 


1 F=.F 8 . 4 , 5 X 4 HGAM=.F 794 ) 

21400000 

001035 

PRINT 120 ,TPRA,EEF 

21500000 

' 0010^5 

120 FORMATI/lX, 22 HRECnVEPY iVi/PuVf ) = /f B / 4 T 7 X 23 Hkl NET I C ENERGY £FF. 

1 =*F 9 « 5 ) 

2160 JO 33 . 
21700000 . 

” 0610 ^ 5 “" 

PRINT 18 *APCA*DELTM,EMAX 1 

21800000 . 

001057 

FORM AT I / 1 X « 22 HC 0 NTRAC T I ON lAtNF/AI == « F 8 « 4 » 27 xiQHnFt TA MAX-.F 8 ^ 4 « 

661657 

36 X 365 

601065 

i 4 X 8 HEP$ MAX=rF 3 . 4 ) 

IFI 06 LTNl,GT.DEiTM|PRINT 119 „ . _ . 

119 FURMAT( 60 X 28 H**«SH 0 CK WAVE IS D£TACHED*« 4 J 
IFIDELTN.GT.DELTNMJGn TO 5 

22 Q 00100 . 
22 iQQQQQ 
2 2 20 3100 
22330000 

301071 

001072 

1 =IF 1 

iF{MC-l) 5 . 20,20 

2 Z 4 C 1090 . 
225 QQQ 0 Q 

001374 

001074 

001376 

20 CONTINUE 



Pl=P 2 

22 6:010 c.. 
22700000 
, 22 . 8 C 0 30 p_ 
22901310 

001077 

Tl=T 2 

OOllOC 

XM=XM 2 N 

230 Q 0 Q 0 Q . 

001162 

_ ti = EliA(KI 

301104 

GO TU 7 



00 U 04 ^ 13 CONTINUE 

OOllOA KM 2 T= 

001112 25 CONTINUE 

oo i ll a XM 2 =( XM 2 T«« 2 «-XM 2 N ^’<=2 )^<'« 5 

001120 AftC~AK«XMNORi^*XM 2 /(XM 2 N*XM) 


00112 A 

001126 

001136 

001141 

001145 


001147 


arca=arca«arc 

ACHG=AC0S( {HC**2fHG#«2-CG**2)/(2.*HC*HG) ) 

AJHG=180.*C-ACHG 

lEIEllAlK KLT*0. ) AJHG=AC‘HG 

AJGH = EP - DELTN „ 

AHJG^iaQ.^C-lAJHG^-AJGH) " 


001153 

001161 

OOllfl 

001 175 

001203 


GJ=HG«(SIN(AJHG)/S IN( AHJG) J 
H j = HO * ( S I N lA J G H > / S T N t AHJG ) ') 

Cj=HC>HJ 

IFieilAlK) .LT.O, )CJ=HJ-HC 
UJ = (GJ*« 2 <-QG=<^« 2)*<'.5 

FPXZP=ACOS{ iQJ*»2VQC*»2-CJ**2) /(2.»aJ»0C) ) 


2330 0 ) 00 . 
234 GJ 00 Q 
23500 090 . 

23600000 
■ZOJQO.ODXL.. 
23B000QQ- 
23900000 
24000000 
2410 3000 
24200000 
24303300 
24400000 
245 30 300 
2460 ppob;_ 
24766000 

24803330 


001214 

001216 


EXZPF=EPX2P 

FXZAVG=0.0 

25000 390 
25100000 

Q01216 


IF (FUA( KM 505, 500.500 

252D30JO 

101221 

505 

AJGHF=EPF 

25300000 

001222 


AHJGF=180-’<'C”U JHG<-AJ'GHF j 

25403000 

001226 


GJF^HG=t'(SIN( ACHG)/SIN( AHJGFM 

25506066 

001234 

001242 


HJF=HG=i= ( SINl AJGHF )/SIN( AHJGFM 
CJF=HJF-HC 

25600 300 
25j000j90_. 

001244 


IF( {EPFH3ELTNKLT,0MCJF = HC-HJF 

25800090 

301251 


OJF=iGJF*^2*-OG**2M^.5 

. 25900000 

>01257 


EXZPF =ACaS((QJF«*2fUC«'*2-’CJF**2)/{2**GJF=«^UC) I 


001270 


IF((EPF<-OELTN) .LF.O.IEXZPF^-EXZPF 

.26100000. 

001274 


EXZAVG' EPXZP-.5«{EPXZP- EXZPF* 

26200300_ 

001300 

500 

continue 

263000QQ 

001300 

JO 6 

FORVIATaOFlO.5) 

26400000 

001300 

001302 


AE0IlP = ATAN2(XM2"N.XM2t) 

26500909 


ABQIIP^^ABQE-AEOUP 

26600000 

001304 


ABEO='ASINiOB*SINIAB’OE*/BE* 

267099QQ 

001312 


AQBE = 180«*C--( ABEO + ABOE ) 

26800000 

001316 


ABI 1P0 = 180.*C'*( AOBEfABOI IP * 

26903030 

001321 


iFUBOriP*3“64.3639 303 

27Q0QQQ0 

001322 

304 

A8I1PQ=a6bE- ABSiABQIlP** 

27100900 

001325 

„_M3 

CONTINUE 

. _ 27200000 

001325 


81 lP=nB»(SiN(A80IlP*/SIN(ABllPQM 

27390 300 

001333 


AABE=ATAN( AE/COI 

27400000 

001337 


T1PI2P== BIIP*CGS<AA 8 E* 

.... 2 7500000 

001342 


AC0I3P = ATAN{T1PI2P/0C* 

27600000 

001346 


AFOC= ATANZlOEtOC* 

2T7C3330 
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001351 

ABUF= ABS(E-AF0C» 

27800000 

001354 

BI 2P= ( B 1 1 P**2^T 1 P 1 2P *«2 1 - 5 

27900 300 

001362 

IF{ DBA0P-l.J200*200,20l 

2800000b 

001365 

001367 

001374 

230 CI 2P = BC -BI2P „ , 

AC0I2P = ATAN(Ci2P/nC> 

GO TO 203 

- uxamisi- 

2820 5000 
28300000 

001374 

001376 

201 Cl 2P = BC tSlZP . 

Ac6l2P=ATAN(ci2P/OCl 

28400330 

28S(lfl0QiL 


203 CONTINUE 


00l‘*03 

EQ^EP/C 

287OQOQ0 

031405 

THETl=^THeTA/C 

28QQOOOQ 

001406 

ueLTNlsOELTN/C _ 

28900000 


_.£PXl=ePX2P/C 

29000000 

001 »11 

EPXIF-EXZPF/C 

29100000 




00141 5 



_001416_ 



001420 

Cai2P=AC0l2P/C 

295000110 

001421 

E0I1P= AEOMP/C 

29600000 

001423 


297QDQQQ 

001 4P 5 





HIBKSEBKRitfrH 

001430 

^80E1==AB0E/G 

iOOOOQQO 

001431 

ABE01=ABE0/C 

30100030 

001433 

DELTN1=0ELTNP/C 

3020.0000 



til !n(i!iltTiW| 

001451 



001453 


3050.0.000 

001454 

RI NFRsR INFR*RR _ 

3Q6QQQ0Q 


PRIMT 114 __ 

307QQQ0Q_ 

00146 2 ~ 
noi 46 P 

114 FORMAT!// ♦ FI nw PRnPFRTTES REFERFNr.FO TO ?sWFPT WFOGF*! 

3a8QOOOQ_ 

PRINT 110 

9nqnnonn 



3100QOOQ 



3110QQQQ 

001466 

PR I NT 111* XM1N*XM1T* X M 2N • XM2T • XM2 * P 1 N FR • R IN FR • T I NFR • GAM 

31203000 

001514 

111 FORMAT! IX. QFin. 51 

31300000^ 

001514 

PRINT 112 _ 

314QOQOQ 

001520 

112 FORMAT(/3X6HDELTAN.5X4HEPSN*6X4HAJOC.5X6HACOI2P^4X6HACOI3P. 

31500000 


14X6HA£OriP.5X5HTH£TA,4X13HOELrAN MAX*3X8HEPSN MAX*3X1HNI 




31700000 

001550 

113 F0RMAT(1X.9F10.5* 171 _ _ 

318Q3QQQ 

001550 

IF(0ELTN1.GT.0ELTM)PPINT 119 

31900001^ 

001556 

PRINT 109, XM,ALP,E1*DEA0R,P1*T1 

32003000 

001576 

109 FDRMATI /IX *3HM1« * F8 . 3* 7X6HSrfEEP= .FB .3* 7X12HR 

. . . _32I003.00. 

LiDGE ANGLE=.F8.3.7X6H0E/A0=.E6.3.7X3HP1=,F8.4,7X3HT1=.F8.2) 32200000 

001576 

PRINT 19» XM2, PIMFR*RINFR*TINFR*GAM 

3230000Q 

001614 

PRINT 120,TPRA.EEF 


001624 

"print 108, ARCA,ePXl,C0r3P 

32500 300 

001636 

13B F0RMAT(/IX*22HC0NTRACTIGN (AINF/AI F8.4t 7X4HEPS^ * F8. 3* 7X11HCRO 

32600000 


ISS FL0W=*F8*3» . . „ 

3270 3000 

001636 

IF(OELTNl,LT.O*)PRINT 508 . EPXl . EPXIF * EXZIAV 

_ _ 32800000 

001651 

538 FOftMAT(/2X.14HEXPANSlCN F AN*. 9X* IDHt EAO WAVE = . F8.4 , 5 XllHFl NAL WAVE 



l = .F8.4t5X9HAVG WAVF='.F8.4) 

333O0IIOO- 

001651 

PRINT 118 __ __ 

33100030 

001655 

118 FORMAT!// * FLOW PROPERTIES RET^ERENCED TO INITIAL XYZ COGROl 

33200000 


INATE SYSTEM* 1 

333C0D00 

001655 

WR = ARCA*SIN< (90.“AL1-COI3PI*C1/COS<"alI*CI 

33400000 

001671 

IF(K.EQ*11PRINT 117»XM2,E11A<K) .EPXl 


001705 

117 FnRMAT(/lX3HM2=,F8.3,3X9H0ELTA XZ= , F8 • 3 , 10X7HEPS XZ-*F8*3» 

33600000 

001705 

IFU^EO^i-ANO.DELfNl.LT.O.lPRINT 5 3 8, EPX 1 ,EPX1F, EXZ 1 AV 

33700000 

001727 

IFlK.eg.lJ PRINT il6,COI3P,WR 

33800000 

001741 

116 FORMAT (/1X28HTOTAL CROSS FLOW lUELTA XY ) = , F8 .3 , lOX 

33900000 


118HGAP RAT IO ( Wl/W2 1“ , F 8 . 3 ) 

. . 34000300 

001741 

IFlMC.eu.ilCO TO 5 

341QQQQ0 

001743 

274 IF(NN1276,606,275 

34200Q00_ 

001745 

606 IF(TYPE(K«'in60l.60l.602 

001750 

276 CONTINUE 

34400000 

001750 

XJOC2=EPXZP 

34500000 

001751 

ACOI2P*COI2P*C 

34600000 

001753 

C0I3P1=C0I3PT*C 

34700000 

001755 

ALP1=I90.-AL1)*C 

34800000 

001757 

FXZP = TANlXJOC2*-ACOI2P3*TAN(ALPl)/(Cr|S(COI3Pl)*(TANiALPll-TAN(C01 

34900000 


i3Pim 

35OQ0J)QO_ 

002000 

EXZP-ATANIEXZPI/G 

35100000 

002003 

XJMUP = EXZPF-‘AC012P 

35200000 

002005 

IF(EXZPF.LT*0*>XJMUP«£XZPF ... ^ 

35300000. 

002010 

EXZPF=TAN( XJMUP 1 *TANI ALPl 1 / (COS (COI 3P11*I TAN( ALPl )-TAN (CO I 3_ 



IPUI) 


3S500000 
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002027 

EXZPF=ATAN(EXZPF1/C 

35600000 

302032 

EXZAVG=exZP-«5=»(FXZP-'EXZPFJ 

35700000 

002036 

0UT2= C AB0E«AE0I1PI/G 

35800000 

002041 

OUTT^COISPTfrOUTZ 

35900000 

002043 

WR = ARCA*SIN( {90«-ALl-GUTT)#CI/C0S<Atl=^CI 

36OOOOO0 

002057 

PRINT 117j XH2,E11A(K),EX2P 

36100000 

002070 

IFiOEiTHl.LT ®0.>PRINT 508,EX2P»EXZPF,EX2AVG 

362d0000 

002103 

PRINT 116, OUTT,WR 

36300000 

002113 

IF(OELTN»GT«DELTNMIGO to 5 

36400000 

002117 

I=I>1 

36500000 

002120 

IFIMC-nS, 5,604 

36600000 

002122 

604 lFlTYPE(K+in601.602,6C2 

36700000 

002125 

275 CONTINUE 

36800300 

002125 

OFT=GOI3PT*C 

36900000 

002127 

XJdC3=EPXZP 

37000000 

002130 

CAPl=190«-ALn*C 

37100000 

002133 

EXZP=TAN{CAPU*TAN {OFTI<'TAN(XJOC3)/ (SIN{OFT I »l TANt CAPll-TAN(DFTI ) ) 

37200 000_ 

002155 

EXZP^ATANC £XZP)/C 

3730Q0Q0_ 

002160 

EXZPF^TAN(CAPU«TAN(OFT)<'TAN( EXZPF )/ (SIN(0FTI *(TAN{CAP1 )-TAN(0FT| 

37400000 


in 

37500000 

002202 

EX2PF = ATANtEXZPF)/C 

37600000 

002205 

EX2AVG= EX2P - ,5*( EXZP-EXZPF ) 

37700000 

002211 

OUTF=COl3PT<-COI3P 

37600000^ 

002213 

WR = ARCA«'SIN({90.-Atl-nuTFnCI/COS(ALl«C) 

37900000 

002230 

PRINT 117, XM2,E11A(K| ,EX2P 

38000000 

002241 

iFIDELTNl.LT.OnPRINT 5 08,EXZP,EXZPF,EXZAVG 

38100300 

002254 

PRINT 116, OUTF,wR 

38200000 

002264 

if(deltn«gt,oeltnm)go to 5 

3830 3300 

002270 

i = H'i 

38400000... 

002271 

1F(MC-‘115.5,605 

38500000 

002273 

605 1F1TYPE<K<-In601,601,602 _ __ __ 

-3.8.6000QO, 

002275 

631 CONTINUF 

38703300 

002275 

NN=~1 _ __ 

38800000 

002276 

CfJI3PT=CQI3PT+C0I3P _ 

38900300 

002300 

AI0I3P^ATAN(TAN(C0I2P*C l*C0S(C0I3PT*Cn/C 

39000000 

002313 

ALP^ATANM TANI ( 90.-ALl-C0I3PT)*Cn/C0S(A1.3l3P4Cn 

_^9100000_ 

002330 

002333 

ALP =90. -ALP/C 

39200300 

K=K+1 .. _ 

39300000 

002335 

FL1=F11ACK-1> 

39400000 



39500Q00_ 

Q 02 36 3 

~ pi = ATAN(TFinC 

39600000 

" 032366 

X.'1 = XM2 

_397000ai]L 

002370 

T1 = T2 

3980QQQ0 


P1 = P2 

39900000 

002373 

0eADR = ABS(TAN(AI0I3P4C)/TANIEl*Cn 

40000300 

002405 

GO TO 7 _ 

40100000 

002405 

632 CONTINUE 

40200300 

002405 

NN=l 

40300000 

302406 

COI3PT = COI3PT4'OUT2 . . 

_4040Q0(1Q_ 

302410 

ALP=ALH'C0I3PT 


002411 

SIGMA=(qO»-ALi)#C 

~ 4C6COOOO 

00241^ 

SIGMAl=C0r3PT*C 

40700000 

002415 

K = K<-1 

40800000 

002417 

E12 = EllAIKI 

40900000 

002420 

E1=ATAN(TAN(E12*C)*(TAN(SIGMA)-TAN(5IGMA1)H'SIN{ $IGMA1)/(TAN{SIGMA 

41000000 


U«tA"N(SIGrtAin) 

41100000 

002446 

El-EI/C 

41200 300 

002450 

XM=XM2 

41300300 

002452 

T1 = T2 

41400000 

002453 

PX = P2 

41500000 

002455 

deadr=o.o 

41600000 

002456 

GO TO 7 

41700000 

002456 

END 

41800000 


SUBROUTINE PEALG (NSTRI 

41900000 


C this SUBROUTINE COMPUTES WAVE PROPERTIES FOR A THERMALLY PERFECT GAS 

^2i^0 3 0 39_. 

300003 

common C,THV,GAMP,RaAH,XMWT,Pl,rL,XM,xMNQRM,TQ 

421000Sl 

000003 

COMMON P2,r2,DeLTN,GAM,XM2N*EP0.RH0R 

4220330Q_ 

303003 

000003 

000003 

common EMAX,DeLTNM,N,KGAS,MC 
COMMON EP0F,TPRA,AR,XMlTtXM2T 
COMMON NDtBUG 

42 300 00C 
424^0 30, 
42500000 _ 

000003 

DIMENSION PC4> ,T14) ,P314) " ' " " ‘ ' 

FGAI T )=l«n GAMP-i.) / (1. 4-<GAMP-1» )♦’{ THEf A/TI**2*£XP( THETA/T) / ( 
2tXP(THETA/TJ-I.i«+2) 

42603033. 
4Z70000g__ 
42 89 JO 33 

00^^ 

FW(T2)=AX'<*2-Q*( A**2-4«*r«T 2)**.5-2.*R*T2-00+GG#R*( T2-TI n-4.*R« 

-^2900000 

000075 

"100076 

000076 

)ooido 

lTHETA»n./(EXP( THETA/T2 )“l. )-ei ) 
EP0F=0» 
ePDAVG=0« 

GC=32«i74l 

GG=4eO#GAMP/ (GAMP-1* ) 

^3003013 
42100 300 
43233003 
42303000 
4 34 CO J03 
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000103 


ft=GC*P8AK/XM^T 

43500000 

'>00106 


RI = P14144«/(RBAR/XMWT=^T1 ) 

43600000 

000112 


THETA^THV 

43700000 

000114 

000115 

XM1=XMW0RM 

ARA^i^O 

43800000 

43900000 

000117 


J*l.. .... - - 

44 000 "30 Q_. 

0QQI2Q 



' 441iHil^L 

000122 

117 

CONTINUE 

44200000 

000122 


DELT==OELTN 

44300 ooe 

000123 


IF(0eLTN*LT,0« SPRINT 118 

44400000 

000132 

118 

FORMAT UX27H»*=S'THERMALtY PERFECT GAS^^^t 14X16HAPPRaX. SOLUTION) 

44500000 

000132 


IF(06LTN.GE«0.)PRINT 6 

4A60000Q 

000141 

6 

F0RMAT(lX27H«««THERMALLy PERFECT 

44700000 

000141 


PRINT 9 

44800)00 

000145 

9 

F0RMAT(38X16HIN FRONT OF WAVE) 

44900000 

000145 


CALL TOTAL I PI , T1 , XM,T 0 , PTl 1 

45000 )00 

000151 


GAMl=FGA(fl) 

45100006 

000154 


XMIN=XM1 

45200000 

000156 


GAM=GAM1 

45300000 

000157 


SOUND =CR*T1*GAM1)44,5 

45400 JOO 

000165 


IFCOELTN.LT.O.) GO TO 446 

45500000 

000167 


EMAX=SQRT( (G AM4-1 .)«( 1 . f (GAM-1 . ) /2. GAM^l. 1/16 •*XMlN4*4 1 ) 

45 60 0000 

000204 


EMAX==(1,/IGAM«XM1N*#2) )*nGAMi>I, )/4.*XMlN4*2-l.i-EMAX) 

45700000 

000217 


EMAX^ASINlSORTIEMAXn 

45800300 

000223 


DFLTNM =(GAM>l*)*XMlN4«2/(2.*(XMiN**2* ( S £N( EMAX l*«2 1-1 . ) I 

45900000 

000234 


0ELTNM=(0ELTNM-1.)*TAN(EMAX) 

46000000 

000241 


DELTNM^ATANa./OELTNM) 

46100000 

000244 


V1-XM1*S0UND 

46200000 

000246 


P I E=3. 141592654 

4630000CL 

000250 


/')2=4 

4640000a 

000251 


IF( UEQ.2)N2=0 

465i)ooaa 

000255 


IF(MC.GT.1)N2=4 

46600000 

000261 


DEL2-DELTN 

46700000 

000262 


I F ( DE L TNM . L T . D E L T N ) D e L 2 ^ 0 E L TN M 

46800000 

000266 


GAN =G AMI 

46900000 

000267 


XM12=XM1*«2 

47000000 

000271 


XM14=XM1**4 

_ 4710000a 

000272 


BB = -(XM12+2* )/XM12“GAN*S1N(0EL2)»’^2 

47200000 

000301 


CC== (2.*XM12 + 1)/XM14>( (GAN<-l.)**2/4,<-(GAN-l.) /XM12)*SIN{DEL2)4*2 

47305000. 

000320 


02=-C0S(DEL2)*^2/XMl4 

47400000 

000324 


EX= (9.*B84CC/2.-BB4*3-27,*D2/2. )/(8B**2-3.*CC)**1.5 

47503300 

000344 


IF(EX.GT.1)EX=1. 

47600000 

000352 


IF(EX.LT*-l)EX=-l 

47700000 

000357 


Ee=ACOS(EX) 

47800000 

000361 


FF = CDS{ (EE«'N2*PlE)/3.) 

47900000 

000370 


G2=(-BB/3,<-( 2^(BB4*2-3.*CC)***5/3. )*FF )^4.5 

48003000 

000406 


EPS==ASIN(G2) 

48106000 

000410 


IF INSTft.EQ. 2. AN0.DEL2 .lt*. ODGO TO 1000 

48203300 

000423 


GO TO lOCl 

48_3000_QO 

000423 

1000 

EPS=90.*C 

48403330 

000425 


GO TO 1004 

40500000 

000426 

1001 

continue 

48603000 

000426 


RR=TAN(EPS)/TAN(EPS-DEL2) 

48700000 

000435 


P«1=1.>GAN» XMi**2*SIN(EPS)**24(l.-l./RR) 

48800000 

000447 

1004 

CONTINUE 

48900000 

000447 


EPl=EPS/C 

49000300 

000451 


M = 1 

49100000 

000452 


N = 1 

49200000 

000453 


NN = 1 

49300000 

000454 


E1^1.7(EXP(THErA/Tl)-l.) 

49403330 

000462 


NT=i 

49500000 

000464 

__ni 

CONTINUE 

49600000 

000464 


KR=-l 

49700000 

000465 


Ul=Vl#SiN(EPS) 

498COOOO 

000470 


A=UUR4T1/Ui 

49900000 

000474 


0=A 

50000300 

000475 


bb=2.64ui=e'*2 

50100000 

000477 


HH=XM1»*2*(SIN(EPS) )**2 

50200000 

000503 


HK=^(2.4GAM1*HH-( GAM1-1. ) ) *( ( GAMl-1 . ) *-2 ♦ 1 / ( ( GAMl U . ) **2#hH ) 

50300000 

000516 


STR=( (-.2/9000. )*T04-1.0 ) #HK 

50400000 

000523 


TM A X= A * * 2 / ( 4 . *R ) - . 0 OObb 1 

50503000 

obosso 


rF(N.FO.l.ANO.NOeBUG.GT.0|PRiNT 223 

50600000 

000544 

\ 22 3 

FORMAT!/) 

50700300 

000544 


IF(NDEBUG.EU.O)GO TO 120 

50800000 

00 05 A 5 


'print 222t tMAXtA,R,UlfTl,VlfEPl,Pkl,PR2 

50 9C 0300 

boost's 

222 

F0RMAT(IX9HSUB REALG,9F12.4) . _ .. .. .. 

SIQQOOao 

000573 

120 

CONTINUE . „ .. — . 

. 5110009.0 

000573 


WMAX::sFW! TMAX) . 

51ZQQQJQ.Q_ 
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300576 

fF(WMAXU40ilOfl43 

, 513Q0C0Q 

000600 

140 


000601 

NT = 2 

51500000 

..M0Q&.Q2 

T21=TMAX 

. 5 160 Q 30 Q 

000603 

T22=T21 

51700000 

000605 

WI-FW(TMAX) 

, 51800300 

000610 

703 T22=^T22-.1«(T21-T1) 

51900000 

00061^ 

W2=PW(T22» 

52000000 

000616 

NT^NT**! 

52100000 

000620 

IF(NT-12I803,8G3»802 

52200000 

000622 

803 lF(W2)703vl9,707 

52300000 

000624 

707 T2i = T22<'.l’MT21-TlJ 

524CJ00C 

000630 

Wl=FWrT21) 

52500000 

000632 

T21=-Wl/(W2-WU»{T22-T21 M-T21 

526 CO 300 

000640 

IF (T21-TMAX) 706*706,708 

527QQQQ0 

000642 

708 T21=TMAX-. 00001 

528C300G 

000644 

706 CONTINUE 

52900000 

000644 

W1=FW(T21) 

53000000 

000646 

T22=T21-.01=«'(T21-T1I 

53100300 

000652 

W2=FW(T22) 

53200000 

000655 

T21=-Wl/(W2-Wll*(T22-r21)fT21 

53300000 

000663 

IF «T21«TMAX> 712*711.711 

53400000 

000665 

711 T21=T 

53500000 

000667 

712 CONTINUE 

53600000 

000667 

W1=FW(T21) 

53700000 

000672 

GO TO 602 

53800000 

000672 

19 CONTINUE 

53900000 

000672 

W1=^W2 

54000000 

000673 

T21=T22 

54100000 

000675 

GO TO 802 

54200000 

000676 

143 CONTINUE 

_ 54300000 

000676 

T21=STR*TH'40. 

54400000 

000701 

IF (T21-TMAX) 300,300, 301 

54500000 _ 

000704 

301 T21=TMAX 

54600000 

000706 

300 T22=T21-10.*XMl 

54700000 

000711 

W1=FW{T21) 

54800000 

000713 

W2=FW(T22> 

54900000 

000716 

T21=-W1/( W2-W1)*(T22-T21J 4-T21 

55000000 

000724 

IF(T21-TMAXm,ll,14 

55100000 

000726 

14 T21=TMAX 

55200000 

000730 

11 T22=T21-2.»XM1 

55300000 

000733 

W1=FW(T21I 

55400000 

000735 

W2=FW<T22) 

55500000 

000740 

T2l=~Wl/tW2”Wl J T22-T21)>T21 . . __ 

5560Q00Q 

000746 

IF (T21-TMAXI141, 141.305 

5570QQ00 

000750 

305 T21-TMAX 

55800000 

000752 

141 T22=T21-.8^XM1 

55900000 

000755 

W1=FW(T21) 

56000000 

000757 

W2=FW(T22I 

56100000 

000762 

T21=-Wl/(W2-W1)#{T22-T21)>T21 

562000X30 

000770 

IF(T21-TMAX>12*12,142 

56300000 

000772 

142 T21^TMAX 

56400000 

000774 

12 T22=:T21-.l*XMl 

56500000 

000777 

W1=FW(T21) 

56600000 

001001 

W2=FW(T22J 

56700000 

001004 

T21^-W1/<W2-WU*(T22-T21)<-T21 

56800000 

001012 

IF (T21~TMAXJ 15,15,10 


001014 

10 T21=XMAX 

57000000 

001016 

15 CONTINUE 

57100000 

001016 

W1-FW(T21I 

57200000 

001021 

802 CONTINUE 

57300000 

301021 

XNT=NT 

57400000 

001022 

XNN^NN 

57500000 

001024 

TRR=T21/n 

57 60 0000 

001026 

T2=T21 

57700000 

001027 

13 CONTINUE 

57800000 

001027 

GAW2=FGA(T2) 

57900000 

001031 

XM2^2.0*T1/IGAM2«T2)*( (GAM1*XMI**2/2.0J <'(GAMP/(GAMP-l.O) T*a.O- 

58000000 


2T2/TU4.THETA/Tl*a./(EXP(THeTA/Tl)-1.0)-l./{EXP(THETA/T2)~1.0in 

58100000 

001067 

IF(XM2i28,28,307 

58203000 

001072 

30'' CONTINUE 

58300000 

001072 

XM2=XM2**,5 

58400000 

001076 

PR2=1«/2«*U 1.I‘GAM2«XM2*«2)-T1/T2*(1«4^GAM1*XM1**2)4‘( ( Ud.+GAMZ* 

58500000 


2XM2*^^2)“T1/T2*( l.t-GAMl + XMl^^a) ) #’4‘2 i <-4. *T1/T2 ) 5> ) 

58600000 

001123 

PR2=l./PK2 

58700000 

001125 

ROR= ( { (GAM2/GAMU*(T2/T1)«( XM2/XM1)«*2-1. » /S IN (EPS ) **2 ) 4- 1 . 

58800000 

001137 

IF (ROR) 10(3,101,102 

58900000 
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001143 

101 CONTINUE 


59000000 

001143 

R 0 R =. 000001 


59100300 

001145 

102 CONTINUE 


592 QQQQ 0 

001145 

R 0 R=R 0 R ***5 


5930000(1 

00115 1 

100 CONTINUE 


5940000 a 

001151 

R 0 R= 1 */R 0 R 


59500000 

001153 

PR 3 = RQR=«'T 2 /T 1 


59600000 

001155 

EP 1 =EPS/C 


59700000 

001157 

IFiNSTkeEQ« 2 *ANO*OEL 2 .CTa* 01 )GO TO 1002 


59800000 

001170 

GO TO 1003 


59900000 

001170 

1002 PK 1 =PP 2 


6000000 a 

001172 

GO TO 25 


60100000 

001172 

1003 CONTINUE 


60203000 

001172 

IF( ABS<PR 2 “PR 1 )”» 00 n 25^25 t 27 


60300000 

001177 

27 N=NM 


60400000 

001201 

NAP=N 


6050000 Q 

001201 

If (M- 2 ) 55 , 56.56 


6060000 Q 

001204 

55 PRY=PR 2 


60700000 

301205 

PRX=PR 1 


60800000 

001207 

IF<NN- 2 I 304 . 303,303 


60900000 

001212 

304 PR 1 =PR 1 >. 5 *(PR 2 -PRI) 


61000000 

001216 

NN = 2 


61100000 

001217 

303 CONTINUE 


612000 QO 

001217 

PRI = PRU-. 01 *(PR 2 -PR 1 ) 


61300000 

001223 

M=MFl 


61400000 

001225 

GO TO 191 


61500000 

001225 

56 S=(PRY-PR 2 >/(PRX-PRU 


61600000 

001231 

PKi = (l«/U«-S)i*(”S*PRX + PRY) 


61700000 

001237 

M = M ~1 


61800300 

001241 

306 CONTINUE 


61900000 

001241 

191 CONTINUE 


620000^ 

001241 

COTD= 1 * 0 /TAN{ OELT) 


62100000 

001244 

6 PS=ATAN(C 0 TD/(GAM 1 *XM 1 * 42 /(PR 1 - 1 . )-l, )) 


6220000 C 

001254 

EP 1 =EPS/C 


62300000 

001256 

If (EPl.LT.l. )GO TO 28 


62400000 

001262 

IF(N- 20 mi. 28,28 


62500000 

001265 

2 B PRINT 112 .EP 1 .T 2 . JU. J.JAA,„NN,NfXM 2 


62600000 

001311 

112 fORMAT(// 2 lH PR DID NOT CONVERGE . 3 X 4 HEP 1 = ,f 10 , 4 , 3 X 3 HT 2 = .F 10 . 4 , 3 X 3 


62700300 


2 HJJ=, I 4 , 3 X 2 HJ=. I 4 . 3 X 4 HJAA=, 3 I 6 .F 10 . 4 ) 


62800000 

001311 

N =100 


62900000 

001312 

Rc TURN 


_ 63000000 

001313 

25 CONTINUE 


63100000 

001313 

P 2 = PRIORI 


A 3200 OOD_ 

001315 

R 2 =R 0 P*K 1 . 


63300000 . _ 

031317 

.MUll-.. 

GO TO 447 - 

. 446 CONTlNUf. . „ 

C C EXPANSION WAVE CALCULATIONS FOR TPG 


- 634 QQ 31 Q 

.^ 63 SQjQLOOO_ 
63 . 600001 } . 

001321 

XNUDF=-DELtN/C 


637 Q 3 Q.Q..Q 


001323 XNUD^XNUDF . _ ^ 

001324 XM|=XMX _ 

00132 6 . T 02 = 10 

0013 2 7 AW = SQRT(<GAM»l , >7 { GAM- 1 . > i 

0013 36 AAy^ = SQRT( (GAM~l, )/( G AM »1 ^ J »(XH I *» 2 ~-l . ) > 

001350 XMU= 90 >~ACQS(i ,/ x'MI i/C 

0 01357 XNUS=AW*ATAN(AAW)/C-(q O«-XMU) _ __ 

Q Q 136 6 _ XN U B = X N U S XN U OF _ 

001367 BW= 1 */AW . 

0013 71 " CW=(GAM>i.r/' 2 V I ' __ 

001 374 _ Dy^=(GAMn»)/( 2 **(GAMrl.J_) 

001377 £>j = -GAM/(GAM-l\) 

001402 XMW=XMI _ . . 

0 01403 _ XNUX^XNUS^'C _ . 

001405 NW = 0 . ~ 


. 6380.0 

639Oiia0O__ 

640 Q_lQilQ 

.6_410100Q 

6420QjQO0_ 

6 44000QO 
6450 3000 
6 460QQ0 0 


. 648 aQQQ.Q._ 

~ 65 0000 Q Q,_ 

65100000 

65201 ) 000 . 


001 4 C 6 444 DO 445 JW= 1,30 


653000 Q 0 


001411 


XNU = AW*ATAN( BW*SgRT(X>^W*XMW’-i. ) J -AT AN ( SQR T ( XMW*XMW-1 * ) ) 

65400030 

001432 


DIFF=XNUX-XNU 

. 65500000 

001433 


If ( ABSIOIf F>->00031)22.22.33 ’ 

656O030C . 

001440 

33 

0M=( XMW/ SQPT(XMW^XMW-i. ) ) * ( 1 * f C W^XMW^X MW f’^'D I Ff 

. 65700000 . 

001452 


XMW=XMW«-DM 

65800 300 

001453 

445 

CONTINUE 

65900000 

001456 


PRINT80 

66000030 

001462 

30 

FOKMAT{30H MAX NO OF [TERATIONS EXCEEQcOJ 

66100000 

001462 

, . . Al 

XMU=ATAN( l ,/SQRT(XMW*XMW-l, n/C 

66200000 

001472 


PR= (i*<-CW*XMW*XMW)**EW 

.. 0630 3000 

001501 


rF(NWl64,64,66 

66400000 

001504 

64 

CONTINUE 

66500 300 

001504 


XMUI = XMU ' ' .. - ~ - 

66 600 00 6 

001505 


prti=pr ■ ■ ^ ' 

66700 300 
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001507 XNUa=XN|iOF 

0Q151O _ _ 1 

001512 xNux= ( xNus ^ fTqat ( nwI^xnud f#e 

001516 7 _ 7 i ^ -X NU8 *C ) 66 ~ 

001522 " 66 CONTINUE “ 

001522 ^ XK2N==XMW 


66B OOOQO 

669000^0 

67000 000 

57100000 . 

672 00000 


001523 

XM2=XM2N 

67400000 

001525 

TR = (T,i-(GAM-I. )/2.*XMlN*«2)/a« <-(GAM-l«)/2.=<^XM2N**21 

675QOOOO 

001537 

T2=T1*TR 

S763000C 

001541 

GAM2-FGAIT2) 

67700000 

001543 

PR1=PR/PRTI 

67800000 

001545 

P2=P1*PR1 

67900000 

001547 

R2=P2*144«/(RBAR/XMWT*T2) 

68000000 

001553 

RH0R=R2/R1 

68100300 

001555 

ePO=XMUl 

68200000 

001556 

EPl-EPD 

68300000 

001557 

EPOF=:XMU 

68400000 

001561 

eP0AVG=EPD-.5^C EPO-EPOF ) 

68500300 

001564 

0=1.0' 

68600000 

001566 

N=JW 

68700000 

001567 

AR=RHOR*(XM2N/XMNORMj*TR**.5 

68800100 

001576 

RNFR=RH0R*XM2N/(XMN0RM*TR^#.26> 

68900000 

001604 

GAMM1=GAM-1. 

69000100 

001606 

TPR = { (GAMn. )*B/{GAMMl*df2. J ) «* ( GAM/GAMMl ( ( G AM^-l . ) / 1 2 . ♦GAM^ 

69100000 


IB-IGAMMU) /GAMMl) 

69200000 

001630 

TPRA=TPRA*TPR 

69300000 

001632 

0ELTN1=DELTN/G 

69400000 


C THIS PACKAGE OVER RIDES PERFECT GAS PRESSURE RATIO AND USES TPG EQ 

69500000 

001634 

GAX = GAMP/<GAMP*-l.l 

69601000 

001636 

ETT=EXPITHETA/T2) 

69700000 

001642 

ETTM1=ETT-1. 

69800000 

001644 

ETTT=EXPITHETA/T02) 

6990J010 

001650 

etttmi=ettt-i. 

70000000 

001652 

ETR=eTTTMl/ETTMl 

70101100 

001654 

ftTR=T2/T02 

70200000 

001656 

THtR=THETA/T2 

70300000 

001660 

THTTR=THETA/T02 

704i):0000 

001661 

EX1=ETT/ETTM1 

70500000 

001663 

EX2=ETTT/FTTTM1 

70601000 

001665 

PT‘PR=ETR«TTTR^*GAX*EXP(THTR*EXl“THTTR*eX2) 

70700000 

001700 

P2=PT1#PTPR 

70800000 

001702 

R2=P2#144./(R8AR/XMwT*T2 » 

70903000 

001706 

RH0R=R2/R1 

71000000 

001710 

447 CONTINUE 

71100100 

001710 

TR=T2/T1 

71 200000 

001712 

XM2T=XM1T«1 IGAM1«T1)/IGAM2^T2) 

71300000 

001717 

XM2F=(XM2T**2i'XM2*#2l«*.^ 

71400000 

00172 5 

PR i NT 7 

71600QQO 

0QI731 

7 FGRMAf(73aXUHBEHlND: WAVEI ^ 

"71603.100. 


_iI 0 I 7 U._ 

■■■OOIJ 35 


001737 

00174_2_ 

0017 ^ 

001754 


call TOTAL tP_2^ T2 iXM^F^_TJ^^ 
..PJR=PT2/_P_TJL_ 


TPRA=TPRA#PTP. 

bELTU=[)ELT/C 

SN»2=(R*T2^GAM2)**.5_ 

AR = R2»XM2»SND2/(R1* X Ml * $ CU N 0 ) 
ARA=AR A»AR . . . 


71 flO O ^00 


0017J5_ 

001 761 

001761 


PRINT 0 
3 FORMAT!// * 
1 EDGE 


fi OW. . PR 0 PT R T I LS ME AS U R E.0 NO R M A L _I_0 . W 6 0 GE LEADING' 


71900000 
72000Q0G ^ 
72 10 OOP G 
72 2QQQ00 
7230.0.QQQ . 
7240000J3:... 
. 7250QOQQ._ 
726QOOOQ_. 


002003 

113 FURMAT(/2X4HMIN=,F6.3,^X7HDELTAN-.F7.3.5X3HP1=*F8-4.6X3HRI=.F8.^, 

IZBOaOM... 

72900100 


26X3HT1=,F8.2,6X4HPT1=,F10.4*6X5HGAM1=,F6.4) 

002003 

PRINT11^,XM2»EP1,P2,R2. T2,PT2.GAM2 

73I0300Q 

73200000 

002025 

114 FORMAT (2X4HM2N=,F6.3* 7X4HEPN= *F7, 3. 5X3HP2= .F8.4.6X3HR?s: . Fft »4 . 


26X3HT2=,F8^2,6X4HPT2=,F10.4,6X5HGAM2=,F6.4) 

002025 

RH0R^P2/P1 

'^3 333030 

002027 

002036 

RNFRVRHbp*XM'2A( XMi#< T2/T1 )**0.26 ) 

PR I NT 1 1 5 , RN FR , AR , P R 1 , R HOR , TR , P TR , N AP 

'73400000 
7? 500000 

002060 

115 FORMAT! lX,5HREYR=,F6.3,5X6HAl/A2=,F7.3,2X6HP2/Pi=.F8.4*3X6HR2/Rl=, 

73600000 


lF9.4,3X6“Ht2/f i=,F8.4,6X4HTPR=,F10.4,9X2HN=,I6) 

73700000 


00206 ^ 

002075 


002075 


IFiOELTD .LT.0.)PRINT 116,___EPPf_FT>i)F,EPDAyG 
116 FGRMAT{/2Xa4HEXPANSlbN F^*,9X,10H WAVE = ,>8.4 ,5 XI INF I NAL WAVE 

. 1 = ,F8,4,5X9H AVG W ~AV E>,H.41 _ 


_A±A±L. 


002077 

0 02100 

002102 

002103 

002105 


XM2N=XM2 
EP0=EP1 
RETURN 
10(9 STOP 
END 


73800000 
73900000 
74 3LO OOOO 
7^V QQ 0Q Q 


T 42 CPOOO 

74300000_ 

7^400000 

_ 7 ^' 50 P^OO 

74600000 


35 



APPENDIX A — Continued 




SUBROUTINE PGAS 

74700000 

c 

THIS SUBROUTINE COMPUTES felAVE PROPERTIES FOR A PERFECT GAS._ 

74800000 

c 

- INPUT= XHlf DELT, N _ _ OUTPUT^ EP 

74900000 

c 

NN= 0 FOR STRONG SHOCK 


c 

NN- 4 FOR WEAK SHOCK 

TOCTI.!>!1M 

000002 


COMMON C,THV*GAMP,RBAR,XMWTfPl*Tl*XM4XMNGRM.T0 


000002 


COMMON P2rT2tQELTN,GAM,XM2NfEPD»RHOR 





75400000 

000002 


COMMON EPDF.TPRA,ARtXMlT,XM2T 


000002 


COMMON NOEBUG 

756Q0QQQ 

000002 




000003 



75800000 

000005 


XM1N=XMN0RM 

75900000 

000006 


EPDF=0. 


OOOOOT 


EPDAVG^O. 

_ 761Q0000 _ 

kiiiitiiiliHHH 


OR = (GAMn.)/{GA«-l.» 

7620D000 

000015 


NN “4 • 

76300000 

000016 


IFIMC.GT.UN=0 

76400000 

000022 


IF (N.EQ«2)NN^0 

76500000 

J00025 


IF(N.EQ.0)NN=4 

76600000 

000027 


IF lOELTNM.LT.OELTN IDELTN=OELTNM 




XM1=XMIN 

76800000 

000033 


13ELT=0ELTN 

76900000 

000035 


PRINT 6 

77000000 

000040 

6 

F0RMAT(1X17H*#*PERFECT GAS***! 

77100000 

0000^0 


PRINT 9 

77200000 

000044 

9 

F0RMAT(3aX16HIN FRONT. OF WAVE) 

77300000 

000044 


CALL TOPGi Pi, TI^XM, TO, PTl) 

77400000 

000050 


IFIDeLTN,.LT.O.I_GO T0 300 

77500000 



P I E=3, 141592654 

77600000 

OOOl53 


XM12=XM1**2 

77700000 

000054 


XM14=XM1*44 


300055 


BS=^(XM12f2. I/XM12-GAM*SIN(DELT)**2 

77900000 



CC=(2.*XM12f l.)/XM 14^1 (GAM M. 1 **2/4. i-I GAM-1 .)/ XM12 )*S INIDELT ) 

78000000 

300102 


DD=-C0STD£LT)#J52/XM14 

78100000 



EX= (9.*BB«CC/2.-B84*3-27.*0D/2,)/(BB*A2-3.*CC1«*1»5 

78200000 

000126 


IF(ex.GT.l)EX-=l. 

78300000 

000133 


IF (EX.LT.-lJEX^-l 


000140 


EE-ACOS(EX) 

78500000 

000142_ 


FF = COSI (EE<-NN*PIE)/3. ) 

78600000 

000151 


GG=(-eB/3.+( 2«( BB4^2-3 .*CC 1 5/ 3. ) *FF ) .5 

78700000 

300167 


EP=ASINIGG) 

78800000 

300171 


FPO=EP/C 

78900000 

00G173 


[F( PPD.GE.89.999)G0 TO 25 

79000000 

300176 


RR=TAN(EP)/TAN(EP-DeLTN) _ 

_ 79100000 

000206 

rii" 

CONTINUE 

79200000 

000206 



79300000 

000220 


PR1^1.«'GAM’«'XM1N**2 4SIN{EP)*^2*( l.-l./RR) 

79400000 

30C 231 


N=NN 

79500000 

000232 


B=XM1N**2*SIN(EP)**2 

_ 796^0000 

303236 


BB = UGAM-1.)*8<-2. >/(2.*GAM*B-(GAM-l«) ) 

79700000 

000244 


XM2N = (B6/(SIN{EP-0£LTNn*42)**,5 

79800000 

000255 


TR = ((2.0*GAM*B-(GAM-1. !*=<'(< GAM-1. 1^8*2* J )/nGA«^l, ) 4*2*8 1 

79900000 _ 

000270 


T2=TR*T1 

80000000 

000272 


P2=PR14P1 

801 00000 

000274 


>R=P2/P1 

00 

o 

o 

o 

o 

o 

o 

000275 


PHOR=RR 

8030 0000 

000276 


EPO=EP/C 

80400000 

000300 


GO TO 502 

80500000 

000301 

500 

CONTINUE 

80600000 

C 


EXPANSION WAVE CALCUtATIGN 

_ 8010030il_ 

000301 


XMI=XM1N 

80800000 

000302 


XNUOF=:-OELTN/C 

80900000 

000304 


XNUO^XNUOF 


000306 


AW = S0RTUGAM«-1.) /{GAM-1 .)) 

81100000 

000314 


AAW = SQRT({GAM-l.)/{GAM4-l.)4(XMI»42-l.n 

81200000 

000326 


XMU=90.-ACOSU*/XMI ) /C 

81300000 

000335 


XNUS=AW*ATAN(AAW)/C-(90.-XMU) 

81400000 

000344 


XNUB=XNUS4-XNU0F 

81500000 

300345 


BW=1./AW 

81600000 

000347 


CW=(GAM-l.)/2. 

81700000 

000352 


0W=(GAM4'1. )/ (2/*(GAM-l . ) ) 

81800000 _ 

000355 


EW=-GAM/IGAM-1.) 

81900000 

000360 


XMW==XMI 

82 00 0 0 00 

003361 


XNUX=^XNUS*C 

82100000 

000363 


NW = 0 

82200000 

000364 


0D10JW=1,30 

8g3Q00QQ 
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000366 


XNU=AW«ATAN( BW* SQRT( XMW*XMW-1 ® ) J 

'-ATAN( SQRT(XMW^XMW-1* I ) 



82400000 

000407 


olfNxnux-xnu 




82500000 

000410 


IF( A8S( 01 FFI- *00001122 r22f 33 




82600000 

000414 

3 3 

OM=( XMW/SgRT(XMW«XMW-l® 1) ♦ ( 1® 4-Cw*XMW^XMW ) *0 1 FF 



82700000 

000426 


XMW^XMW^'OM 




82800000 

000427 

10 

continue 




82900000 

000432 


PRINT80 




83000000 

000435 

80 

FJRHATOOH MAX NO OF ITERATIONS 

EXCEEOEOJ 



83100000 

000435 

22 

XMU=ATANC 1«/SURT(XMW*XMW-1. ) ) /C 




83200000 

000445 

000454 


PK= (l.4-CW«XMW*XMW)**EW 
IF«NW^55*55.66 





833001^0. 

834aQaao_ 

000456 

,55 

CONTINUE 




83 5000:^0 

000456 


XMU1=XMU 




8 3601 >00 

000457 


PRT1=PP 




837000QO 

000461 

000462 

00046'4 


XNUO=XNUDF 

NW=NW4-1 

XNUX= ( XNUS 4-FLOAT { NW LAXNJJO ) 

•• 



838>0l0i3_. 

83900000 

8^006000 

000470 

000474 

66 

IF (XNUX-XNUB«C)444*444,66 
CONTINUE 

-• -- 



a^ioom, 

8^200003 


00047A 

000475 

60647 7 

0 00512" 

^0051^ 

000516 


XM2N=XMW _ 

PR1 = PR/PRT1 _ _ .. I- J 

TR = U * f( G A M~ 1 * ) / 2 . * 1 N « <1 » » ( G AM- / 2 . * X« 2N* # 2 » 

T2 = Tl:^Tfi ■ J ‘ ^ ^ 

P2 = P1«PR1 

RHQR=PR1/TR " ‘ 



000517 
„ 000521 
"00052 2 
A® 526 
0005 30 
000532 


EPO-XMUl 

5POF=XMU 

EPD AVG^EPO-. 5=* ( EPO-EPOf ) 

8=1.0 

N=JW 

50 2_€ 0 N T I N UH 


0^00532 

000 541 

000547 


000551 

IMSUL 


0006 00 

" 000^6 

00061^ 

000612 

00061 6 

000620 


AR=RHOR<^( XM2N/XMNORM)«TR=^'>i=.5 

RNFR = RHGR«XM2N/(XMNOR'm*TR**.26) _ _ _ _ 

GAMMl = GAM-l . ” J'-J - ^ 1 IT " 

T P R^ ( ( G AM 4- 1 . I * 8 / ( G AMM1*6 + 2 . ') 7* * ( G A M/G A MMi ) * f( G AM >1 V ) / ( 2 • *G AM * 
lfl-<GAMMl)))«#(l./GAMMl) __ _ 

2 T= X_M ij^ ( Tl/T2)»*.^ ' 


84^0 0 O OP 

35oa6opo_. 

35100660 

85200000 

85300000 

85466000 


0006 21 

0QQ6 22 

■)00624 

000624 

J0Q6 30 


XM2F^(XM2T«4 2 4-XM2N**2 l«+.5 
PRINT 7 

7 FORMAT ( /33 X 1 IHBEH IND_ WAV E » 

CALL T0PG(P2,T2jXM 2F ,T02t_PT2 > 

tpka=tpra<^tpr 

ue LTN1= DELTN/C 


. acr=o. 

GAM 1.= GAM 

GAM2=GAM 

Rl=Pl*144,/(R8AR*U/XMWTr 
R2=P2*14A./(RBARxtT2/XMWT )_ 


85500000 
85600 06 c 
85700666 
85.80 OOPP_. 

85960060 
86000000 
86100000 
86200000 
86 30J_OOa 
86400000 
865PPPP0. 

" 86JG mo . 

868Q0QQQ 
. 86 9C0300 

mmop. 

8Ja6gopq_ 


000636 

PRINT 8 

8730QQQQ 

.000^4 2 

8 FORMAT!// # FLOW PROPERTIES MEASURED NORMAL TO WEDG>= LEADING 

8740 30-30 


1 EDGE ♦) 

87500300 

000642 

PRINT 3.XMN0RM,0£LTN1.P1.RI.T1.PT1.GAM1 

87600000 

000664 

3 FORMAT(/2X4HMlN=.F6.3.4X7H0ELTAN=.F7.3,5X3HPr^,J=8.4.6X3HRl = .F8_.4^ 

R770000Q 


26X3HT1=.F8.2«6X4HPT1 = .F10.4.6X5HGAM1 = .F6’.4I 

87800000 

000664 

PRINT 4, XM2N.EPD,P2.R2.T2.PT2*GAM2 

87900000 

000706 

4 F0RMAT(2X4HM2N=,F6.3*7X4HEPN=.F7.3,5X3HP2=.F&.4,6X3HR2=:.F3.4. 

88000030 


26X3HT2=*F8«2»6X4HPT2=»F10*4*6X5JH_GAM2=.F6«4> 

38100000 

000706 

PRINT ?.RNFR.AR.PR1.RHGR.TR,TPR,N 

88200009 

000730 

2 FQRMATI1X,5HREYR=.F6.3,5X6HA1/A2=«F7.3*2X6HP2/P1=.F8.4.3X6HR2/Ri=. 

88300000 


lF8.4.3X6HT2/Tl=»F8.^t6X4HTPR=,FlQ«4.9X2HN=.I61 

8840 0000 

000730 

IF(OELTN1.LT.O.|PRINT 5t EPO^EPOFtEPOAVG 

38500300 

000743 

5 F0RMAT{/2X,14HEXPANSI0N F AN*. 9X . lOHLEAO WAVE =. F8.4 , 5X11HF I NAL WAVE 

88600000 


l=.F8.4.5XqHAVG WAVE=.F8.4)- 

88700 a 30Q 

000743 

RE TURN 

88800000 

000744 

END 

88900000 


SUBROUTINE TOT AU P »T , YM , T Tqt7 W 

89000000 



THIS SUBROUTINE COMPUTES TOTAL PROPERTIES FOR A THERMALLY PERFECT G.4S 

_ 891QQQQQ 

boooio" 

COMMON C TTHVyGAMP*R6AR* XMWT.PI .Tl.XM.XMNORM.TO 

89200000 

000010 

COMMON P2tT2«0ELTNtGAM* XM2N*EP0 sRH 0R 

'8930Q0QQ 

000010 

COMMON EMAX»nELTNM,NjKGASf MC 

89400300 

ooooio 

COMMON EPDF,tPRA,AR,XM*lT,XM2T 

89500000 

000010 

COMMON NDEBUG 

39600000 

000010 

FGACT r=1.4-(GAMP-i. )/(1.4-(GAM“P-la I«(THETA/T)*«2*EXP(fHETA/t) /{ 
2EXPCTHETA/T)-1.)*’4'2) 

.a^mQQoo... 

89800000 

000044 

FXMCTOI = SQRTUA*TO-T,/GAM>«'B’^(l./(EXP(THV/TO)-l.V-En 

89900000 
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APPENDIX A - Continued 


000072 

00007^ 

000077 

000077 

THeTA=THV 

PRINT2 \ ^ 

2 FORMAT! 49X, ♦TOTAL PPOPERT I J 
PKINT3 

aQ.hOQSlQ.Q— 
90200000 
9030 0000 

^ooyo3 

3 FORMAT(8XlHM,9X6HMCH£CK,7X2HTOf 9X4HT/ r0,eX2HP0,l0X4HP/P0» 5X1HK, 
IUX4HR/R0*8X3HGAM,6X1HN » 

90400000 

90500000 

000103 

M=1 

906JOOOO 

000104 

GAM=FGA(T) 

90700000 

000 111 

RHa=P4144./{R8AR«T/XMWT> 

90BCQ030 

0001 i 5 

7*/{GAM«T) 

90900 000 

000120 

2««THV/(GAM*T) 

91000000 

000122 

ET= exP(THV/T) 

. 9UCOOCO 

000130 

E=1»/(ET-1,) 

91200000 

000133 

TOl * T*M,#-CGAM-l.)*yM*^2/2.) 

91300000 

000140 

T02=T01-50, 

91400000 

000142 

8 XM1=FXM(T01) 

91500000 

000144 

Xv^2= FXM(T02I 

91600300 

000147 

TO=(yM-XMl)/(XM2-XMU*(T02-foiH-T01 

91700000 

000156 

XM3=FXM(T0) 

91800000 

000161 

IF( A8S(YM-XM3H -0001) 10.10 » 5 

91900000 

000166 

5 T31=TC 

92000000 

000167 

TOZ^TO-IO. 

92100000 

ooom 

N=Nfl 

92200000 

000173 

IF(N-25)9,9,10 ... 

923D0300 

000175 

9 GO TO 8 

92400000 

0C0176 

10 CONTINUE 

92500000 

000176 

ETO=EXP(THV/TO) 

926COOOO 

000202 

ETOl=ETO-l- 

92700000 

000204 

ET1=ET-1. 

_ 92800000 

000206 

APR = EXP(tTHV*6T/CT4ETl ) )-(THV4ET0/{T04ET0i n ) 

92900000 

000225 

PR- ET01*(T/T0)*^{GAMP/(GAMP--1. ))*APR/ET1 

.. 93Q.OOOOO 

000236 

RHOR = ET01^IT/T0>**I1-/*GAMP-1.))4APR/ET1 

93100000 

000247 

TR-T/TO 

93200000 

000251 

PO-P/PR 

933QQ0QQ... 

000253 

TTOT-TO 

934QQQQ0 . 

000253 

PRINT1,YM,XM3*T0,TR.P0.PR.RH0.RH0R.GAM.N 

93500000 

J003D3 

1 FORMAT(lX,2fl2.5.F12.3,F10.5,F13.3,F13.9,F12.6.F13.9.F9.5.I5» 

93600000 

000303 

RETURN 

93700000 

000304 

END 

93800000 _ 


SUBROUTINE TOPG ( P , T , YN , TTOT , PO 

. 939Q000Q 

c 

THIS SURftnijflNF "c'gMPUTFS TOTAL PROPERTIES FOR A PERFECT GAS . 

94000000 

000010 

COMMON C.THV.GAMP.RBAR.XMWT.PI .Tl.XM.XMNORM.TO 

941QQ0QQ - 

000010 

common P2.T2, DELTN. GAM. XM2N.EPD.RHGR 

942QQQQQ..,-. 

000010 

COMMON HMAX.OELTNM.N.XGAS.MC 

943QQ0QQ . 

000010 

COMMON EP0F,TPKA.AR,XM1T,XM2T 

944QQQ0Q 

000010 

COMMON NOEBUG .. . 



PRINT2 

94600000 

2 FORMAT! 49X.4T0TAL PROPFRTIES*» 

947GOOOO 

”000013 

PRINT3 

948QQQQQ-. 

000017 

3 FORMAT! 8X1HM,9X6HMCHECK. 7X2HT0.9X4HT/ TO. 8X2HP0.10X4HP/P0.9X1HR. 

949QQQQQ 


1 10X4HR/R0,aX3HGAM.6XlHN J 

95000000 

000017 

RBAR-1545.31 

95IQi)QQQ_ 

000020 

XMwf=28.964A 

95200000 

000022 

RH0-P4i44./(RBAR4T/XMWT) 

953Q0QQQ 

000030 

A - (GAM-i.l/2.*YM*42 

95400000 

000034 

TR - l./(l.FA) 

95500000 

000036 

PR = <1./(1.4-A))««(GAM/(GAM-1.)) 

95600000 .. 

000045 

RHOR = {l./(l.«-Ai)4*a./(GAM-l.n 

9570000Q _ 

000055 

TO-T/TR 

95800000 _ 

000057 

PQ=P/PR 

959QQQQQ 

000060 

TTOT=TO 

96000000 

000061 

XM3=YM 

96100000 

000062 

PRINT1.YM,XM3.T0.TR,P0.PR.RH0.RH0R.GAM 

962QQQ0Q 

000110 

1 FQRMAT!1X,2F12.5,F12-3.F10*5.F13.3.F13.9,F12.6.F13.9.F9.5.I5) 

96300QQQ 

000110 

RETURN 

964iI0000 

loom 

END 

_ 9650.0000 _ 
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APPENDIX A - Continued 


**** SAMPLE CASES 



CASE NO. 1 — FOR AN IN iTtAl SWEEP ANGLE OF 48 DEGREES. THIS CASE COMPUTES 
FLOW PROPERTLES ACROSS 3 SUCCE SSt VE SHOC K WAVES. THE WAVES 
ALTERNATE TYPE A^ TYPE B/ AND TYPE A. WIT H THE FLOW TURNING 
ACROSS^CM SHOCK W~a Ve~ EQUAL TO 6 DEGRE ES WHEN MEASURED IN 


I NPUT FOR CASE " j 

^ 48. 6,0 , 02200 

S 6 . 6 . _ ,, _ _ .... ■ ■ . ■ 

, . -■ , ... ■ 


„ CASE. NO, 2 — TJilS IS ALSO A REJECT hkS CASE^ WITH AN INi T lAL S WEE P ANGLE 

OF 48 DEGREES. THE f[rST1 V;AvL ti A S HOC K FOLLOWED BY AN 

. . . . EXPANS I ON FAN BOTH. W AVES ARE TYPE A. ANO T HE FLOW TURN ING 

- - . ACROSS EACH WAVE IS. E.QUAL_TQ. MEASU RED ^JQjE 

XZ_ilLAN£^__ . . 


l-NPUT FOR CASE NO. 2 _ ^ 

2 ^ 6.0 48. 6.0 0. _ U4 _ _ 02200 

6^ - 6 . _ _ _ 1 _ 

-i.. ...... .1. : _r:. . :i.. :: 


CASE NO. 3 THIS CASE IS A REPEAT OF CASE NO . I, BUT FOR A" THERMALLY 
PERFECT GAS AND AT A FLIGHT DYNAMTC PRESSURE OF” 1000 PSF7 

INPUT FOR CASE NO. 3 

3, 6, P „48s, 6, 0_ 0^ . 2775 402. 22200 

6, G. 6. 

1 . - 1 , 1 . 
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APPENDIX A - Continued 
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APPENDIX A - Continued 
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APPENDIX A - Continued 


42 


TOTAL CROSS FLOW (DELTA XV)^ 5^084 GAP 



CASL-N0__,2 WAVE NO 


APPENDIX A - Continued 
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»*^PERF£CT GAS*»« 


APPENDIX A - Continued 


31 O Cvl 

I— •4' 
O N. ivj 

Z ♦- r-l 
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TOTAL CROSS FLOW (DELTA XYI= ^159 GAP RATIO (Wl/W2)s *973 



»»*thermally perfect ga$*»* 


APPENDIX A - Continued 
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TOTAL CROSS FLOW (DELTA XY)- 2>9Q7 GAP RATIO (m/W21? 2,654 




APPENDIX B 


INLET PERFORMANCE FOR LOW REYNOLDS NUMBER 
AND UNSTARTED OPERATION 

In order to evaluate the inlets operatir^ characteristics better, one test was made 
with the Reynolds number reduced by a factor of 1/3 to 3.3 X 10® per meter. The result- 
ing static -pressure distributions within the inlet are found in figure 61 where they are 
compared with the high-pressure tests. Also given in figure 61 are the results of a test 
in which the model failed to start because of blockage produced by a survey rake. 

The low Reynolds number tests had only a small effect on inlet operation as deter- 
mined by static -pressvire level chaises produced by slight alterations in shock-wave 
positions. In the side passage these effects were seen near the cowl in figures 61(f) 
and 61(h), whereas the center -passage effects are observed in figure 61(k). The high- 
pressure level indicated on the foreplate (fig. 61(a)) is in error because of the diffi- 
culty in measuring the extremely low pressures. 

Alterations in static-pressure levels for the unstarted choked inlet were observed 
throughout the model. High pressures were found on the foreplate (fig. 61(a)), top sur- 
face (figs. 61(b) and 61(c)), sidewalls in front of the struts (fig. 61(d)), and on the cowl 
(figs. 61(f) and 61(g)). The pressure level in the center passage was low and uniform 
and indicated the absence of shock waves (figs. 61(i) and 61(k)). 
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APPENDIX C 


CAPTURE MEASUREMENT DATA 

The primary conclusions from the capture measurement data (Mach number and 
capture parameter contour maps) taken downstream from the struts are presented in the 
report in figures 53 and 54. The purpose of this appendix is to include the data from 
pitot and static -pressure surveys used to produce those results. Like the data taken at 
the inlet throat, the program used to analyze these data and generate the contour maps 
is discussed in appendix B. 

The static -pressure distribution around the walls of the inlet at the capture meas- 
urement station is given in figure 62 where X is the peripheral distance around the area 
as defined in the sketch. The static probe Survey data are found in figure 63 which also 
includes the wall values of f^ure 62. A nominal value of 80 percent was selected for the 
total-pressure recovery limit at the capture measurement location downstream of the 
struts. When the static surveys were combined with the pitot surveys of figure 64, the 
computed total-pressure recovery was above the imposed limit of 80 percent only in the 
very small region indicated by the dashed line in figure 63(d). The data are seen in con- 
tour map form in figures 65 and 66, and a total-pressure recovery map is given in 
figure 67. 
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TABLE I.- THEORETICAL FLOW FIELD PROPERTIES USED FOR FIGURE 5 
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TABLE n.- THEORETICAL END EFFECTS USED FOR FIGURE 5 


(a) Top surface (b) Cowl ^6jjy = 0°) 


Bay 

M 

P/Pl 

PtA,i 

^pitot/^1 


Bay 

M 



Ppitot/Pl 

2 

4.92 

3.07 

0.939 

97.1 


2 

4,99 

2.82 

0.942 

91.7 

3 

4.46 

5.22 

.910 

136.0 


3 

4.28 

6.48 

.907 

156.0 

4 

4.13 

7.79 

.888 

175.0 


4 

3.67 

14.40 

.877 

256.0 

5 

3.92 

10.25 

.883 

208.0 


5 

3.35 

22.30 

.867 

333.0 

6 

3.77 

12.53 

.879 

235.0 


6 

3.04 

34.90 

.854 

432.0 

7 

4.67 

4.26 

.933 

122.0 


7 

4.48 

5.17 

.930 

136.0 

8 

4.25 

6.80 

.903 

161.0 


8 

3.86 

11.30 

.894 

222.0 

9 

3.98 

9.53 

.891 

199.0 


9 

3.31 

23.60 

.870 

344.0 

10 

3.72 

13.58 

.879 

248.0 


10 

2.88 

43.40 

.846 

484.0 

11 

3.56 

16.47 

.871 

277.0 


11 

2.62 

63.60 

.829 

592.0 

12 

4.46 

5.22 

.910 

136.0 

I 

12 

4.28 

64.80 

.907 

156.0 

13 

4.13 

7.79 

.888 

175.0 

1 

13 

3.67 

14.40 

.877 

257.0 

14 

3.92 

10.27 

.886 

208.0 


14 

3.35 

22.40 

.870 

334.0 

On plow* 

5.33 

1.98 

.968 

73.3 







*In front of bay number 2, 
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TABLE m.- CORRECTED THEORETICAL FLOW FIELD PROPERTIES 


USED FOR FIGURE 27 


Bay 

M 

P/Pl 

Pt/Pi 

^xz 

®xy 

A/Ai 

PpitotAi 

1 

6.00 

1.00 

1.000 

0 

0 

1.00 

46.8 

2 

5.07 

2.54 

.922 

6.83 

2.00 

1.80 

85.2 

3 

4.46 

5.12 

.891 

0 

3.18 

2.96 

133.0 

4 

3.97 

9.44 

.870 

6.00 

5.40 

4.44 

196.0 

5 

3.70 

13.61 

.866 

2.00 

6.52 

5.66 

246.0 

6 

3.45 

19.14 

.862 

6.00 

8.27 

7.07 

302.0 

7 

4.66 

4.09 

.912 

10.83 

3.19 

2.56 

116.0 

8 

4.07 

8.37 

.879 

4.00 

5.20 

4.11 

182.0 

9 

3.59 

15.55 

.858 

10.83 

8.14 

6.16 

265.0 

10 

1 

11 

12 

3.24 

25.52 

.842 

4.83 

10.28 

8.45 

357.0 

4.38 

5.59 

.879 

0 

3.44 

3.13 

141.0 

13 

3.90 

10.19 

.859 

-6.00 

5.71 

4.65 

204.0 

14 

3.64 

14.62 

.857 

-2.00 

6.87 

5.91 

256.0 

*15 

3.40 

20.40 

.852 

-6.00 

8.66 

7.34 

313.0 

16 

3.24 

25.28 

.835 

4.83 

11.30 

8.38 

354.0 


*For large center strut. 
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TABLE IV.- CORRECTED THEORETICAL END EFFECTS USED FOR FIGURE 27 
(a) Top siirface =* 4°) (b) Cowl a 0°) 


Bay 

M 

p/Pi 

PtA,! 

Ppltot/Pl 

2 

4.86 

3.25 

0.921 

100.0 

3 

4.39 

5.61 

.891 

142.0 

4 

4.07 

8.26 

.870 

180.0 

5 

3.88 

10.60 

.866 

210.0 

6 

3.73 

12.90 

.862 

522.0 

7 

4.58 

4.46 

.912 

123.0 

8 

4.17 

7.34 

.879 

168.0 

9 

3.87 

10.60 

.858 

209.0 

10 

11 

12 

3.63 

14.60 

.842 

254.0 

4.33 

5.95 

.878 

146.0 

13 

4.03 

8.56 

.859 

183.0 

14 

3.84 

11.10 

.857 

216.0 

*15 

3.70 

13.40 

.852 

242.0 

16 

3.69 

13.80 

.835 

248.0 


Bay 

M 

p/Pi 


Ppitot/^1 

2 

4.86 

3.25 

0.921 

100.0 

3 

4.18 

7.22 

.996 

166.0 

4 

3.57 

15.80 

.858 

267.0 

5 

3.28 

24.30 

.849 

348.0 

6 

2.97 

37.80 

.835 

447.0 

7 

4.35 

6.52 

.902 

162.0 

8 

3.68 

13.40 

.867 

240.0 

9 

3.10 

31.10 

.829 

399.0 

10 

11 

12 

2.69 

56.30 

.801 

551.0 

4.11 

8.08 

.873 

180.0 

13 

3.50 

17.40 

.845 

283.0 

14 

3.22 

26.60 

.838 

368.0 

*15 

2.90 

41.10 

.822 

465.0 

16 

^ 

2.63 

59.90 

.783 

562.0 


*For large center strut. *For large center strut. 
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Figure 1.- Airframe-scramjet-engine integration. 
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Figure 2.- Inner module of scram jet ei^ine concept (ref. 8). 




Vehicle uaderhody 



Figure 4.- Inlet model sketch and Mach 6 shock diagram. 



Flow- bay - (see table I) 










Center strut 



(b) Throat relief next to cowl. 




(e) Cowl leadii^-edge design. 
Figure 6.- Concluded. 
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Figure 7.- Off-design shock wave systems (xz-plane). Sweep = 48°. 
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Figure 8.- Model photograph. 



L-72-7494 


(b) Front view. 
Figure 8.- Concluded. 
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Left sidewall 



F^ure 9.- Model schematics. All dimensions are in centimeters 










Figure 10.- Side and center strut dimensions as measured in the xz-plane. All dimensions are in centimeters 




Cowl leading-edge positions 



Stations are measured relative to sidewall leading ec^e; cowl leading edge shown at 
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Figure 13.- Static orifice locations. H » 19.05 cm. 





Sidewall, side^passage throat 
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Figure 13.- Continued. 
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(d) Cowl. 

Figure 13.- Continued. 
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Figure 13.- Concluded. 
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Figure 14.- Thermocouple locations. H * 19.05 cm; 24-gage iron-constantan. 
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Figure 20.- Inlet model mounted in test facility 
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Disturbanc 
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(a) Three-strut configuration (Capture = 81 percent) 


FLOW 


Wm 


(b) No center strut (Captiu:e ® 92 percent) 
Figure 24.- Oil-streak photograph. 














Vehicle imderbody 



SECTION A-A (x2-plane) 

Figiire 26.- Corrected Mach 6 shock diagram. 










Figure 27.- Corrected shock -wave detail in vicinity of struts. 
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Throat 



s7H 

Figure 30.- Static -pressure distribution. Sidewall; Y/H * 0.43. 
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Throat 

















s'7H 


Figure 35.- Static -pressure distribution center passage. Side strut; Y/H = 0.43 





s"/H 

Figure 37.- Static -pressure distribution center passage. Side strut; Y/H = 0.88. 



Cowl lea^iing edge 



c/H 

Figure 38.- Slatic-pressure distribution. Cowl (center line). 
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Q Taken for initial inlet configuration 

and shifted to match new strut position 
□ Modified configuration 
Large center strut 



Figure 39.- Static -pressure distribution. Cowl (Z -direction). 
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Figure 40.- Static -pressure distribution. Throat (center -passage side strut). 
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Figure 41.- Static-pressure distribution. Throat (center-passage center strut). 
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Figure 42.- Oil -flow study 














Figure 45.- Throat pitot pressure surveys. The solid symbols represent 
theoretical boimdary -layer calculations. 





Figure 45.- Concluded, 
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Figure 46.- Throat static -pressure surveys. The solid symbols represent 
data faired to the wall static pressure. 
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Mach number 



(b) Center passag( 
Figure 48.- Mach number 
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(b) Y/H = 0.28. 
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Figxire 63.- Static -pressure surveys at capture station. Mj^ = 6.0, 
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Figure 63.- Concluded. 
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Figure 64.- Concluded. 
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